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Abstract 

Based on the analysis of 6.031xl0 7 /T ± — > TtPifin* decays, the NA48/2 collaboration has recently 
determined the S-wave nn scattering lengths oo - o-i with high precision. In addition, the scatter- 
ing length «2 has been independently measured, although less precisely so. The present article 
discusses in detail one of the theoretical frameworks used in the data analysis. 
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1. Introduction 

Half a century ago, Budini and Fonda [1] investigated the threshold singularities in A' ± — > 
n c decays and showed that nn rescattering generates a cusp in the partial decay rate dF/dE„± . 
The strength of the cusp is determined by the amplitude for the reaction n + n~ — > n°n°. Budini 
and Fonda provided an analytic formula for the cusp behavior and pointed out that, in principle, 
this decay allows one to measure the nn S-wave scattering length aq - a 2 , where the a; denote 
the scattering lengths of definite isospin / = 0, 2. There were only a handful of K — > 3n decays 
available in those days, and it was impossible for the authors to actually determine ao, a 2 in 
this manner. The method was then forgotten and rediscovered 45 years later by Cabibbo in his 
seminal work on the interpretation of the cusp detected in data on K — > 3n decays, collected 
by the NA48/2 collaboration B ■ 

In the last decade, spectacular progress has been achieved concerning the knowledge of nn 
interactions, in theory as well as in experiment. As for theory, the scattering lengths were pre- 
dicted with percent level accuracy 04j,|5|], 

a = 0.220 + 0.005, a 2 = -0.0444 ±0.0010, Oq - a 2 = 0.265 ± 0.004 , (1.1) 

within a framework which combines Roy equations |@] and chiral perturbation theory (0, 01. 
On the experimental side, progress was mainly achieved through the analysis of three specific 
processes. First, pionium decays into two neutral pions allow one to measure \o.q - a 2 \ fl fioll . 
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Experimental results have been published by the DIRAC collaboration lUlll . Second, using Wat- 
son's theorem and numerical solutions of the Roy equations lfl2l[l3ll . it is possible to determine 
flo and fl2 from Km decays. Experiments with spectacular statistics have been carried out in this 
channel also in the last decade lfl4l [l5 , 161. Last, the above mentioned cusp in 



decays has been investigated through the analysis of 6.031xl0 7 such events, and a high preci- 
sion value for ao — a 2 is now available lfl7ll . It confirms the chiral prediction very nicely. Indeed, 
combining K e \ decays and the result from the cusp analysis, the most recent publication by the 
NA48/2 collaboration quotes 

ao = 0.2210 ± 0.0047 stat ± 0.0040 syst , 
a 2 = -0.0429 ± 0.0044 stat ± 0.0028 syst , 
ao - a 2 = 0.2639 + 0.0020 stat ± 0.0015 syst . (1.2) 

It seems fair to say that the precise values Eq. (11.2b could only be obtained through a combined 
effort of experiment and theory. Indeed, relating experimental data to the scattering lengths is 
nontrivial affair [18, Ijl 2(J, and a precise theoretical description of the cusp behavior in the 



cl 



amplitude for A' ± — > n "n n* turns out to be quite difficult 



In their most recent data analysis [17], the NA48/2 collaboration makes use of decay ampli- 
tudes constructed along two different frameworks. In the first one ^ |2j]], the structure of the 
singularity at the cusp is investigated using unitarity, analyticity, and the cluster decomposition. 
In addition, an approximation scheme is used, which consists in expanding the decay amplitude 
in powers of nn scattering lengths. The latest work [21] retains effects up to order (scattering 
lengths) 2 and omits radiative corrections. 



The main purpose of the present article is a detailed description of the second method [22, 



23j,|24J]. It uses a Lagrangian framework, which automatically satisfies unitarity and analyticity 



constraints and allows one to include electromagnetic contributions in a standard manner 02411 . 
In order to retain the possibility of an expansion in powers of scattering lengths, which is a very 
convenient concept, a non-relativistic framework is invoked that has already proven to be useful 



in the description of bound states [25], see Ref. [26] for a review. The formalism has recentl 



also found applications in various, mostly cusp-related studies 12711 . such as rj — > 3n 128. 
and 77' — > TJ7T7T 130 [ decays, as well as in near-threshold pion photo- and electroproduction on 
the nucleon 113 ll I32I1 . The amplitudes for K — > 3n decays have been evaluated within other 
frameworks as well. We shall come back to these in Sect. [9] 

The two main difficulties we were faced with to evaluate the pertinent amplitudes are the 
following. 

• The huge statistics available require a very precise theoretical description of the decay 
amplitude, which allows to determine the scattering lengths from data with high precision 
as well. In particular, it is mandatory to include the effects of real and virtual photons 
["radiative corrections"]. 

• The decay amplitudes can be calculated only within a certain approximation. It is therefore 
mandatory to set up a power counting scheme that permits to quantify the neglected effects 
in an algebraic manner. 

As far as we can see, quantum field theory is the only method that allows to satisfy these two 
requirements. Independently of the method used, one is in addition faced with the problem that 
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the decay amplitudes for K — > 3n are beset with leading Landau singularities [33, 341 135L 1361137 , 



decay 
39]Q 



This does not come as a surprise, because we are dealing with unstable particles here. 



Our article is organized as follows. In Sect. |2] a covariant, non-relativistic quantum field the- 
ory framework is set up and applied in Sect. |3]to nn scattering, which plays a crucial role in the 
analysis of final-state interactions in K — > 3n decays. In Sect. [4] we present the Lagrangian for 
K — » 3n decays and determine the tree contributions, while one-loop effects (two-loop effects) 
are treated in Sect. |5](Sect. |6). In the latter section we also evaluate the pertinent two-loop inte- 
grals in the non-relativistic theory in some detail. Vertices with six-pion couplings are shown to 
be of no relevance at the present accuracy in Sect. Q while we present the structure of the com- 
plete amplitude at two-loop order in Sect. [8] A comparison with other approaches to K — > 2>n de- 
cays is provided in Sect. [9] and a summary and concluding remarks are given in Sect. [10] We have 
relegated technical aspects of our work to various Appendices. In particular, in Appendix A we 



provide the complete results for the two-loop amplitude. The explicit expressions for the perti- 
nent two-loop integrals are given in Appendix B and their holomorphic properties are discussed 
In order to make certain that we were not lead astray with the non-relativistic 



Appendix C 



framework, we have evaluated several of the emerging two-loop graphs in standard relativistic 
quantum field theory as well. In particular, in Appendix D the case of relativistic two-loop in- 
tegrals with equal pion masses is considered, while loops containing different pion masses are 
investigated in | Appendix E| Landau singularities are consi dered in | Appendix F| and a compar- 



ison with the non-relativistic loop integrals is performed in Appendix G Finally, we comment 



on the decomposition of the amplitudes in singular and non-singular parts, as used in Ref. H21H . 
in Appendix H 



2. Covariant non-relativistic framework 

In this section, we construct a non-relativistic effective field theory for the two-particle sec- 
tor, which will later be used to describe pion-pion scattering. In order to simplify the discussion, 
here we first give a formulation for the case of a single self-interacting scalar field cD(x) with 
mass M. In the following, the theory will be extended to describe three-particle decays. The 
main difference to the conventional non-relativistic formalism is that our approach as developed 
here — despite the fact that we deal with a non-relativistic framework — yields explicitly in- 
variant two-body scattering amplitudes^ The present form turns out to be very convenient to 
describe three-body decays, since, in particular, the location of the singularities in the two-body 
subsystems (which are generally not in the rest frame) coincide with the exact relativistic values 
to all orders in the non-relativistic expansion. In the conventional setting, this can be achieved 
only perturbatively, order by order. Such a course of action is not wrong, but very cumbersome 
and certainly not elegant. 

In a non-relativistic theory, one assumes that the momenta of all particles are much smaller 
than their masses. The question of the domain of applicability of the non-relativistic description 
is very subtle and depends on the dynamics of the particles considered!! This issue will be 



'in the following, we adhere to the notation used in Refs. 137114011 and use "anomalous threshold' as a synonym for 
"leading Landau singularity". 

2 For an elementary introduction to the essentials of the conventional setting in the context of the nn scattering, see, 
e.g., Ref. [ 26] and references therein. 

3 We remind the reader that we are designing an approach with the intention to deal with the rescattering of pions in 
the final state of K — » 3;r decays. However, as it is easy to check, the maximal momenta of the pions in this decay are of 
the order of the pion mass and thus the applicability of the non-relativistic theory in this case is not a priori clear. 
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discussed in detail in the following sections. For the moment we merely assume that using the 
non-relativistic framework can be justified, and proceed with its precise formulation. 

Constructing an effective field theory implies setting certain counting rules. We formalize 
the definition of the non-relativistic domain by introducing a generic small parameter e and 
postulating the following counting rules in this parameter. 

• The mass M is to be counted as 0(1); 

• all 3-momenta are counted as |p ; | = 0(e); 

• consequently, the kinetic energies are counted as Ti = p® - M = 0(e 2 ). 

The modified approach differs from the conventional one in two aspects. First, the usual 
non-relativistic propagator has a pole at p° — M + p 2 /2M that corresponds to the non-relativistic 
dispersion law. The higher-order corrections to this dispersion law are treated perturbatively, so 
that the results of loop calculations with the non-relativistic propagators can be made Lorentz- 
invariant only approximately, at a given order in the expansion in M~ l . In the modified frame- 
work, these higher-order corrections to the one-particle propagator are summed up, leading to 
the relativistic dispersion law p° = w(p) = -y/p 2 + M 2 . 

The second modification is related to the matching of the non-relativistic and relativistic 
theories. In the conventional non-relativistic approach, the matching condition between two- 
body scattering amplitudes, which fixes the values of the non-relativistic couplings in terms of 
the parameters of the underlying relativistic theory, reads (see, e.g., Ref. [26]) 

4 

|~[(2w(p;)) 1/2 r AIS (p3,p4;pi,p2) = 7s(p 3 ,p 4 ;pi,p 2 ), (2.1) 

i=i 

where the subscripts NR (R) stand for the non-relativistic (relativistic) amplitudes. The addi- 
tional factors (2w(p,)) 1/2 for each external leg account for the different normalization of the 
non-relativistic and relativistic states. 

Since Eq. ( 12. U is not explicitly covariant, the matching condition is different in different 
reference frames. For simplicity, we consider matching at threshold, where the relativistic am- 
plitude is merely a constant 7«(p3, P4; pi , P2)| thl . = A in all reference frames. On the other hand, 
the non-relativistic amplitude at threshold in different frames, obtained by expanding Eq. (I2.lt 
in momenta, is given by 



7V«(P3,P4;P1,P2) 



thl - 4M 2 16M 4 



P 2 + . . . , (2.2) 



where P denotes the total 3-momentum of the particles 1 and 2 at threshold. Let us now limit our- 
selves to the tree approximation and suppose that one wishes to write down the non-relativistic 
effective Lagrangian that reproduces Eq. (12.2) in this approximation. Such a Lagrangian would 
have to consist of an infinite tower of operators, whose couplings are determined by a single con- 
stant A. This again means that Lorentz invariance can be taken into account only perturbatively, 
order by order in the expansion in M~ l . Once more this procedure, albeit formally correct, looks 
rather awkward and renders higher-order calculations cumbersome. 

In the particular case considered above, it is clear that the problem disappears if we arrange 
the non-relativistic theory in a way such that the overall non-invariant factor on the left-hand 
side of Eq. (12. U disappears. This can be achieved by a non-local rescaling of the non-relativistic 
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field <D(x) — > V2W O(x), where W = VM 2 - A . In the modified theory, the normalization of 
the 1-particle states is given by the relativistic expression (p|q) = 2w(p)(2n) 3 6 3 (p - q), whereas 
the matrix element of the free field operator between the vacuum and the one-particle state is 
normalized to unity: (0|<t>(0)|p) = 1. The Lagrangian with the rescaled field is given by 

L = &2W(id, - W)0> + C(O i ) 2 <D 2 + . . . , (2.3) 

where C denotes the lowest-order (0(1)) four-particle coupling and the ellipsis stands for four- 
point interactions with derivatives (these terms count as 0(e 2 ) and higher). Note that the above 
Lagrangian is non-local since it contains square roots of a differential operator. Heavy-particle 



number should be conserved in the non-relativistic theory by construction (see, e.g., Ref. [26]). 
The matching condition in tree approximation yields T^r = 4C = A in all reference frames or, in 
other words, the truncation of the higher-derivative terms in the Lagrangian does not render the 
non-relativistic amplitude non-invariant. 

The propagator in the rescaled theory is given by 

i{0\T®(xW(y)\0) = — L (2.4) 

J (2nf 2w(p)(w(p) - po - le) 

From now on, we shall not display ie in the propagators explicitly. Our next goal will be to 
demonstrate that 

• in the above effective theory, the counting rules in e, established at tree level, are not 
destroyed by loop corrections; 

• the amplitude is explicitly Lorentz-invariant also in the presence of loops. 

Since all loop diagrams in the non-relativistic approach to the two-particle sector can be ex- 
pressed through the elementary bubble integral 



J( P° P) = | — (7 5) 

J (2n) D i 2w(I)2w(P - l)(w(l) - f)(w(P - 1) - P° + Z°) ' v ' ; 

we concentrate on this integral below. The quantity J(P°, P) is calculated in an arbitrary reference 
frame characterized by the total 4-momentum P = (P°, P). Calculations are done in dimensional 
regularization, D is the number of space-time dimensions and d — D - 1 . 

To begin with, we note that naive power counting predicts J(P , P) to be of 0(e d ^ 2 ): each 
propagator in the integral scales as 0(e 2 ), while the integration measure, with d momenta and 
one energy integration, contributes a power e d+2 . It can be checked by straightforward calculation 
that the loop given by Eq. ( 12.5b violates this power counting prediction. This is a well-known 
phenomenon, caused by the presence of the heavy mass scale M in the integrand. In order to 
circumvent this problem, usual Feynman rules in the effective theory have to be supplemented 
by some additional prescription that annihilates the high-energy contribution (coming from the 
integration momenta of order M) to the Feynman integral. We choose a particular prescription 



r d D i 



1 



referred to as the "threshold expansion" 114111 (see also Refs. [26, 42] for a focused discussion of 
the issue within non-relativistic effective theories). The prescription defines how the square roots 
present in the particle propagators are handled in the calculations. According to the prescription, 
one expands the integrand in a Feynman integral in inverse powers of the heavy scale M, inte- 
grates the resulting series term by term in dimensional regularization, and finally sums up the 
results. 
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Let us see how this prescription can be realized in practice. We perform the integration over 
Z° in Eq. (EJ) by using Cauchy's theorem and rewrite J(P°, P) as 



j ^-J£? ' - - 



i i 

4w(l)w(P - 1) w(l) + w(P - 1) - P° 



Further, performing the shift 1 — » 1 + P/2, we transform the integrand by using the identity 

1 1 1 ) 



J_f L 

V a Wh I P° - w„ 



4\V a Wb I P° -W a -Wb P° + W a + Wb P° + W a - Wb P° - W a + Wb J 

1 1 



2P0 12 _ (lp/po )2 _ q 2 ' 



(2.7) 



where w fl = w(P/2 + 1), w b = w(P/2 - 1), q\ = s/4 - M 2 , and s = P 2 = (P ) 2 - P 2 . 

Next, we investigate how the threshold expansion affects the result of the integration. Ac- 
cording to power counting, w a — M, Wb — M = 0{e 2 ) and P° - 2M = 0{e 2 ). Consequently, 
expanding the last three terms in Eq. (12. 71 in powers of e generates polynomials in the integra- 
tion momentum 1. Recalling that integrals containing only powers of 1 vanish in dimensional 
regularization, we finally conclude that using the threshold expansion enables us to rewrite the 
integral in Eq. ( 12.6b as 



J(P°, P) J ( ^ )d 2pQ p _ _ . (2.8) 



(1P/P°) 2 -^ 2 

It is seen that the heavy scale M has disappeared as a result of the threshold expansion. 

In order to perform the integral in Eq. (12. 8) . we choose the first axis along the momentum 
P, so that P = (|P|, 0) and 1 = Rescaling l\ — > l\ P / t/s and doing the momentum 

integration, we finally arrive at 

J(P°, P) = J(s) = — [1 - — j + 0(d - 3) . (2.9) 

Note also that J(s) coincides with the imaginary part of the relativistic one-loop integral. 

The non-relativistic amplitude in the absence of derivative couplings is given by the bubble 
sum 

r WR (p 3 ,p 4 ;pi,p2) = 4c + 8C 2 y(^ + i6C 3 y(i) 2 + ... . (2.10) 

The inclusion of derivative couplings is straightforward. We restrict ourselves to order e 2 , where 
the real part of the relativistic scattering amplitude of two identical particles can be parameterized 
in terms of two constants A and B, which are related to the S-wave scattering length and the 
effective range, 

Re r s (p 3 , p 4 ; Pi , p 2 ) = A + B(s - AM 2 ) + 0(e 4 ) 

= A + B(p lP2 + P3P4 - 2M 2 ) + 0{e A ) . (2.1 1) 

The term with two derivatives in the non-relativistic effective Lagrangian can be directly read off 
from Eq. (12.1 11 1, leading to 

L {2) = D{(W&W® f ® 2 + VO 1 VO 1 "® 2 - M 2 (<D' ) 2 <D 2 ) + h.c] , (2.12) 
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with AD - B + C 3 /(64tt 2 M 2 ). A gain, this Lagrangian is non-local, and the factor W should be 
expanded in actual calculations. 

We summarize our main findings up to this point. 

1 . If the threshold expansion is applied, the one-loop integral J(s) counts as 0(e), so the loops 
do not violate power counting. 

2. The insertion of derivative couplings does not violate power counting. 

3. J(s) is explicitly Lorentz-invariant. Since the tree-level amplitude is also invariant, so is 
the scattering amplitude at any given order in the perturbative expansion. This statement 
is trivial for non-derivative couplings only; one can easily ensure that it still holds in the 
presence of derivative couplings in the scattering amplitude. 

4. According to Eq. d2.91 >. J(s = AM 2 ) = 0. This means in particular that the coupling con- 
stant C is proportional to the scattering length to all orders in perturbation theory. In 
general, the coupling constants of the non-relativistic theory are expressed through the 
effective-range expansion parameters in the two-particle sector. If in the following this 
non-relativistic Lagrangian is used to evaluate pion-pion rescattering effects in the three- 
particle decay in perturbation theory, the result will be written in terms of these parameters. 
This property constitutes the major advantage of the non-relativistic approach as compared 
to a relativistic framework. For example, the decay amplitude calculated in chiral pertur- 
bation theory is given as an expansion in the quark masses, not in the scattering lengths; if 
one attempts to extract the values of the scattering lengths from the data on three-particle 
decays, the accuracy of the former representation may not suffice. 

The generalization of this method to higher orders and to the case of non-identical particles with 
different masses can be performed in a straightforward manner. 

We conclude this section by a remark about the terminology used. We still refer to the above 
framework as "non-relativistic," albeit the energies and momenta of particles in this approach 
obey relativistic dispersion laws. In our naming scheme, "non-relativistic theory" denotes a 
theory in which explicit antiparticle degrees of freedom are absent (respectively, are included 
in the couplings of the effective Lagrangian), and the number of particles is preserved in each 
interaction vertex. 

3. Non-relativistic approach to mt scattering 

3.1. Lagrangian and scattering amplitude 

Now we apply the modified non-relativistic framework to the nn scattering amplitudes. The 
masses of charged and neutral pions M„+ = M„ and M n o are taken to be different, but virtual 
photons are not included at this stage. Due to the inelastic coupling of the 7r + 7r~ and the 
channels, a consistent power counting requires the quantity = M 2 - M 2 to be counted as 
0(e 2 ). We consider the following five physical channels in n a n b — > if-rf 1 : (ab; cd) = (1) (00; 00), 
(2) (+0; +0), (3) (+-; 00), (4) (+-; +-), (5) (++; ++). The Lagrangian takes the form 

5 

£™ = Yj °1 2w ±0' 5 < " w ±)®± + *o 2W b(»i - Wb)O + J] JU , (3.1) 
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where O, is the non-relativistic pion field operator and W± - y/M 2 . - A, Wo = y Af 2 - A, with 
A the Laplacian. We furthermore introduce the notation 

(®n)» = {Pn)^n , = (PnMPn)v®n , {V n \ = (W„, -/V) , 

= (P^<Dj , (^V = (^(PDv®! ' = «V) , (3.2) 

for n = a, b, c, d, in order to write the interaction Lagrangian in the form 
JU = XiQ^l^b + h.c) 

+ - (.®tr®l]k®a) fl ® b - + /;.c.} + . . . , 0.3) 

with /z, = - 1 (M 2 + M 2 , + M 2 + M 2 ,) , where 5, denotes the physical threshold in the z'th channel. 
Explicitly, h\ = 2M 2 , h 2 = 2M n M JI o, h 3 = 3M 2 - M 2 , h A = h s = 2M 2 . The ellipsis stands 
for terms of order e 4 . The low-energy constants C,, D,-, E, are matched to the physical scattering 
lengths below. To simplify the resulting expressions, we have furthermore introduced the scaling 
x\ - Xf, - 1/4, X2 — X3 — x\ — 1, u\ — 113 — us — 0, U2 — U4 = 1. Finally, note that we do 
not discuss local six-pion couplings here, which can potentially play a role in the rescattering of 
three-pion final states; we will comment on these in some detail in Sect. [7] 

Evaluating the non-relativistic scattering amplitude up to order e 2 with the use of the above 
Lagrangian, we obtain 

TZ = 2Coo + 2Dqq(s - soo) + 2C 2 00 Joo(s) + 4C 2 x J + -(s) 

+ 2C 3 m (Joa(s)f + SC oC 2 x J 00 (s)J + -(s) + 8C^C + _(/ + _(i)) 2 + . . . , 

r+° = 2C +( , + 2D +0 (s - 5 +0 ) + E +0 (t - u) + 4C 2 +Q J +Q (s) + 8C 3 +Q (J +l) (s)) 2 + ..., 

T* R = 2C, + 2D x (s - s x ) + 4C x C + -J + -(s) + 2C x C 0Q J 0Q (s) + 8C 2 + _C x (J + -(s)) 2 
+ 4(C + -C,C 00 + Cl)J + .(s)Joo(s) + 2C 2 lQ CAJoo(s)) 2 + ..., 

Tj£ = 2C+_ + 2D + _(s - 5 + _) + E + .{t - u) + 4C\_J + -(s) + 2C 2 x Jqo(s) 
+ 8C*_(/ + _(*)) 2 + &C + -C 2 J + 4s)Joo(s) + 2C 2 x C 00 (J 00 (s)) 2 + ..., 

T£ = 2C ++ + 2D ++ (s - 5 ++ ) + 2C 2 ++ J ++ {s) + 2C 3 ++ (J ++ (s)) 2 + ..., (3.4) 

where 5, denotes the threshold in the pertinent channel Soo = 4M 2 , s + o = (M„ + M„o) 2 , s x = 
s + - = s ++ = 4M 2 . Further, in order to make the expressions more transparent, we have modified 
the notation according to 

{C U C 2 ,C 3 ,C 4 ,C 5 } -» {Coo,C +0 ,C v ,C + -,C ++ } , (3.5) 

and similarly for D, and Finally, J a b(s) denotes the generalization of the loop function J(s) 
given by Eqs. i2.5\ and ( 12. 9t to the case of unequal masses M a and Mb, which is obtained from 
Eq. (O by replacing w(l) -> w a {\), w(P - 1) -» w b {V - 1) and WajbQs) = (M 2 l b + k 2 ) 1/2 . This 
function is equal to 

i , 4o 2 ,(s) X(s,M 2 M 2 ) 

Jab(s) = — v ah {s) , v h (s) = = 5 , (3.6) 

\6n s s z 
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with the usual Kallen function A(a, b, c) — a 2 + b 2 + c 2 - 2{ab + ac + be). 

For the purely elastic channels (+0) and (++), Eq. (13.41 > can be directly compared to the 
effective-range expansion of the relativistic tttt scattering amplitudes 



T l R = 32tt^(2Z + l)?;P/(z) , i = +0, ++ 

1=0 

Re 4 = Ai + B iq 2 + 0(e 4 ) , Re t\ = q 2 A^ + 0(e 4 ) , (3.7) 



where the t, are the partial waves of the respective (physical) channel i, A,-, B,-, A^ the corre- 
sponding S-wave scattering length and effective range as well as the P-wave scattering length, 
and z - cos 8 is the cosine of the center-of-mass scattering angle. In the isospin symmetric limit 
at M„ = 139.6 MeV, which we denote by C, etc., the matching relations are of the form 

/_ C 3 \ C 3 

2C + n = C++ — \6na2 , 2 [D+o — - = D++ — — - = Anbo , E+a — \2na\ , 

I (32ttM^) 2 / 02nM n ) 2 

(3.8) 

with the / = 1 P-wave scattering length a x = (0.379 ± 0.005) x 10 _1 Mr 2 in addition to the S- 
wave scattering lengths of definite isospin / = 0, 2 quoted in Eq. (11. U l4j,|5|]. For the remaining 
channels, the situation is more complicated due to the coupling of the n + n" and channels 
with their different thresholds. In the isospin limit, with both thresholds coinciding, one has 
matching relations also for these similar to Eq. (13.8) . and the coefficients C, and £, are related to 
the nn scattering lengths according to 

\6n - 167r - 8tt 

Coo = — (flo + 2a 2 ) , C x = — (a 2 -a ) , C+_ = — (2a + a 2 ) , E+- - \2na\ . (3.9) 



The relations for the D, are similar to those displayed in Eq. (13.8) . However, taking the pion 
masses to be different, the threshold behavior for tttt channels of total charge zero is rather of the 



form (cf. Ref. lUJ) 

32tt f = 0? + iv+-0 + /voor' + v+-Voo$ (3.10) 

(for the case of an S-wave), where a', . . . , 5' are real and analytic functions in the low-energy 
region on the real axis, < s < 16M 2 . The above parameterization is an analog of the effective- 
range expansion for the tttt amplitudes in the presence of multiple thresholds. Although we do 
not make use of these expressions in the following, for illustration we show the corresponding 
functions for the (00) channel, up to the order displayed in Eq. d3.4l >. i.e. 0{a 3 e 2 ), which read 

a 00 = 2(c 00 + 4D o4) " 7J^( C o>oo + 4C + _C%_)) , 

o00 1 n 2 00 1 n 2 r00 ' n r 2 n1n 

P =Tn C *' 7 = 8^ Co0 ' 6 = ~32^ C00C - (3 - U) 

We note in passing that recovering the isospin limit from these coefficients is a delicate procedure: 
it is obvious that the second term in the g 2 -expansion of a 00 will not coincide with the effective 
range in the isospin limit (up to normalization) due to the presence of the term <5 00 . The ft', y', 
6' are not independent, though, but also these can be expressed in terms of the parameters of the 
effective range expansion. The constraints that enable us to determine these parameters can also 
be derived from (multi-channel) unitarity, which must be obeyed by the tttt scattering amplitudes. 
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The effective field-theoretical approach that we use automatically incorporates unitarity order by 
order. 

In order to take isospin breaking in the leading nn effective range parameters into account, 
we consider the effective 0(p 2 ) Lagrangian of chiral perturbation theory 

F 2 

Xeff = —{dpU&'tf + X U Jf +x'U) + e 2 C(QUQU*) , (3.12) 

with x — 2B diag(m„, mi) , Q = diag(2, — 1)/3 , and U being the standard 2x2 unitary matrix of 
the pion field. Further, at this order, 

M 2 , = M 2 = B(m u + mi) , M 2 K = M 2 + . (3.13) 

Note that evaluating isospin breaking only at 0(p 2 ) is consistent with neglecting virtual photons 
in the non-relativistic framework. In order to carry out the matching of the relativistic and non- 
relativistic theories at higher orders, it is necessary to include photons. At leading order in chiral 
perturbation theory one finds 14311 

C o, + o, ++ = Coo, + o, ++ (l-77), C X = C,( l+7?/3), C + _ = C+_(l + tj) , (3.14) 

where rj = A„/M 2 = 6.5 x 10~ 2 . Equation ( 13.141 ) shows that the threshold amplitudes, which 
occur in the scattering length expansion, are affected by substantial isospin breaking corrections. 
Isospin-breaking corrections in the effective ranges and P-wave scattering lengths (or, more pre- 
cisely, in the non-relativistic couplings D, and £,) do not contribute at the accuracy considered 
here. They have been investigated in Ref. 02911 and were found to be tiny. 

Generalizations to higher loop orders in nn scattering are straightforward. However, they 
are irrelevant for the calculation of K — > 3>n decays to two loops, as performed below; in fact, 
the corrections of 0(C y ) in the matching relations for the effective ranges Eq. (13 .8b are already 
beyond the order needed for our purposes, as they constitute 0(a 3 ) effects and only enter K — > 2>n 
at three loops. 

3.2. Convergence of the non-relativistic expansion 

We wish to briefly discuss the convergence of the non-relativistic representation of the nn 
amplitude as given in Eq. ( 13.4t . As an example, we show real and imaginary part of the 7 = 
S-wave amplitude in Fig. Q] over the kinematic range accessible in K — > 3tt decays, i.e. for 
4M 2 < s < (Mk — My} 2 . We compare to the Roy equation solution lfl2ll matched to chiral 
perturbation theory 0|. We observe that at such low energies, already the tree amplitude 
in Eq. ( 13.41 i. consisting merely of scattering length and effective range term, gives a very good 
description of the real part, see Fig.[T](left); including the two-loop corrections hardly changes the 
amplitude at all. This is easily understood by simplifying the 1 = S-wave two-loop amplitude 
with the matching relations analogous to Eq. (13.8b according to 

Re t° = flo + b l)q 2 + J^f/ 2 ^ , (3.15) 
M 2 (M 2 + q 2 ) 

so it differs from the tree amplitude only by terms of 0(a 3 q 4 ) and higher. A further improvement 
requires the introduction of a shape parameter oc coq 4 . On the other hand, the description of the 
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Figure 1: Convergence of the t = 0, / = tttt partial wave in the non-relativistic expansion; the left panel shows the 
real, the right panel the imaginary part. In both cases, the full lines represent the Roy equation solution 1 12] matched to 
chiral perturbation theory [4, 5]. In the left panel, the dashed line is the non-relativistic tree-level amplitude, while the 
dash-dotted curve represents the two-loop representation. In the right panel, the dashed curve is the one-loop amplitude 
according to Eq. j3.41 . while the dash-dotted is the improved one-loop representation, see main text. 

imaginary part, given purely by the one-loop contributions in Eq. ( 13.4ft . misses the phenomeno- 
logical amplitude completely, see Fig. Q] (right). This is also easy to explain: even in the sense 
of (only) perturbative unitarity, the imaginary part of 0(a 2 e) corresponds to a real part in the 
scattering length approximation. We have to improve the imaginary part to 0{a 2 e 3 ) and partial 
0(a 2 e 5 ) by the replacement C, — > C, + Z),(s - *;) in the one-loop contributions of Eq. (13.4b ; this 
improved one-loop representation is then a very good approximation to the true imaginary part 
at the energies relevant for K —> 3n, see also Fig. Q] (right) 

Finally, we wish to anticipate the use of the tttt amplitudes discussed above in the cusp anal- 
ysis of K — > 3n. The decay rates will be expressed in terms of the couplings C,, D,-, . . . that 
are related to the physical nn scattering amplitudes. Once these quantities have been determined 
from data (in practice, one may decide to use some of the parameters, for instance effective 
ranges or P-waves, as input, employing their theoretically predicted values), they can be related 
to the S-wave scattering lengths using the corrections displayed in Eq. ( 13.14b . With radiative 
corrections applied [24], these relations must be adapted accordingly. 

4. Non-relativistic approach to K + -* 3n decays 

4.1. Kinematics 

In this and in the following sections, we develop in detail a non-relativistic framework 
for K — > 3n decays by means of the particular channels K + (Pk) — > Ji i) {pi)n () {p2)7t + {pi) and 
K + (Pk) — > Jt + {p\)n + {p2)7t~{p?,). The method can straightforwardly be extended to neutral kaon 
decays - for early applications of the method, see Refs. 11221 l23l 12411 . 

The kinematical variables are defined as follows, 



Si = {P K -Pif, M 2 =p 2 , 
3s° ±M 2 K + Ml+2M 2 , 3sl±M 2 K + 3Ml, T^rf-Mi, j = 1,2,3, (4.1) 
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with 



Si + s 2 + s 3 =M\ + J]m 2 . (4.2) 



Further, assuming that i, j, k are all different, we set 

a2 ^ A(M 2 K ,Mf,(Mj + M k ) 2 ) 
\M\ 

In the center-of-mass frame Pk = (Mk, 0) , 



Qi = Py + P*, Q-=P°j + P° k - (4.3) 



0= M 2 K + Mf- Si A(M 2 K ,M 2 , Si ) 
Pi 2M K ' P '' 4M 2 K ' ; 

As usual, the label i — 3 is always assigned to the "odd" particle, i.e. to the n + in the neutral 
mode K + — > 7r°7r () 7r + and to the ii~ in the charged mode K + — > 7r + 7r + 7r~. Note that the values of 
the masses M-, are channel dependent. 



4.2. Lagrangian and tree amplitude 

A non-relativistic approach to describe decays K — > 37r can be justified if the typical kinetic 
energies T, of the decay products are much smaller than the masses. This can be achieved by 
considering a world where the strange quark mass is taken to be smaller than its actual value. 
Then, a consistent counting scheme arises if one introduces a formal parameter e (the same as in 
the two-particle case) and counts T, as a term of order e 2 , the pion momenta as order e, whereas 
the pion and kaon masses are counted as 0(1). From Pi = Mk - one concludes that 

the difference Mk - 2, Mi is then a quantity of order e 2 as well. In addition, as mentioned in the 
previous section, the pion mass difference A n must also be counted as 0(e 2 ). The effective field 
theory framework, which we construct below, enables us to obtain a systematic expansion of the 
amplitudes in e. For sufficiently small m s , the expansion in e is expected to work very well. 

Together with e, our theory has another expansion parameter, namely a characteristic size of 
the nn threshold parameters, which we denote genetically by a. In particular, the amplitudes in 
the non-relativistic framework are given in form of an expansion in several low-energy couplings 
Q, £>,, Ej, which can be expressed in terms of the threshold parameters of the relativistic nn 
scattering amplitude. We expect the expansion in a to converge rapidly because of the smallness 
of the scattering lengths. These two expansions are correlated: because one-loop integrals are of 
order e, adding a pion loop generated by a four-pion vertex increases both the order in a and in e 
by one. A consistent power counting is achieved: to a given order in a and in e, a well-defined 
finite number of diagrams contributes. 

Increasing m s to its physical value again, convergence in the e-expansion is not a priori evi- 
dent, because TJMi can become as large as 0.4, and the corresponding maximal momentum |p;| 
is then not much smaller than the pion mass. However, let us note that the non-relativistic frame- 
work is only used to correctly reproduce the non-analytic behavior of the decay amplitudes in the 
kinematical variables s\, S2, S3, and to thus provide a parameterization consistent with unitarity 
and analyticity - a trivial polynomial part in the amplitudes can be removed by a redefinition 
of the couplings in the Lagrangian. In addition, from the analysis of the experimental data one 
knows [21] that in the whole physical region the real part of the decay amplitude can be well 
approximated by a polynomial in si, S2, S3 with a maximum degree 2. We interpret this fact as 
an experimental indication for a good convergence of the e-expansion for the quantities we are 
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interested in. In the following, when we also study the behavior of the decay amplitude in the 
complex energy plane away from the real axis, we therefore understand the low-energy region 
(or, equivalently, the non-relativistic region) to be determined as a strip enclosing the real axis 
from si = (Mj + M^) 2 to s, = {Mk - Mi) 2 , and going slightly beyond the boundaries. The width 
of the strip should be smaller than the hard scale set by the pion mass squared. 

We now proceed with the construction of the Lagrangian framework. Aside from the La- 
grangian SL„ n displayed in Eq. ( 13.11 ) that describes nn final state interactions, we need the La- 
grangian £,k which generates genuine K — > 3n decays, such that the complete Lagrangian is 

-C = -C K +-Cnn- (4.5) 

At order (scattering lengths) 2 , the amplitudes are given by 

[K + -> ttW] , 
[K + -> ttVtT] , (4.6) 

with obvious notation. The Lagrangian X.k is now chosen such that the tree-level amplitudes up 
to and including 0(e 4 ) become (cf. Ref. U^, Eqs. (4.6), (4.7)) 

Mf e {s u s 2 , S 3 ) = X Q + Xi(sa - s°) + X 2 (s 3 - s° ) 2 + X 3 ( S1 - s 2 ) 2 , 

M^ e (s u s 2 , s 3 ) = Y + Yi(33 - sj) + T 2 (s 3 - s + ) 2 + Y 3 ( Sl - s 2 ) 2 . (4.7) 

We assume T-invariance and a hermitian Xx, as a result of which the couplings Xj, T,- are real. 
Expressing s* through p9 with the use of Eq. 04.41 1. these expressions are equivalent to 

Af(ii,s 2 , *3) = Go + Gdp° 3 - M K ) + G 2 (p° 3 - M n f + G 3 {p° l - p\f , 

Mf e ( Sl , s 2 , S3) = H + H x (pi - M„) + H 2 (p° 3 - M n ) 2 + H 3 {p\ - p° 2 f , (4.8) 

where 

Go = X + ((Mjc - M K ) 2 - s°) X x + ({M K - M„) 2 - s°) 2 X 2 , 

Gi = -2M K Xi - AM K ({M K - M K ) 2 - s°) X 2 , G 2 = AM 2 K X 2 , G 3 = 4M|X 3 , 

H = T„ + ({M K - M n ) 2 - sj) Tj + ((M* - M„) 2 - s+f T 2 , 

Hi = -2M^Tj - 4M* ((Mjc - M,) 2 - s + ) Y 2 , H 2 = 4M|T 2 , ff 3 = 4M 2 T 3 . (4.9) 
From the expressions Eq. (14.8) one may read off the pertinent Lagrangian, 

£ K = tf2W K {id, - W K )K + ^ + h.c) + ^ (tf{W + - M,)0 + <D 2 + h.c) 

+ y (^ f (W + - M^) 2 + Og + /lc.) + G 3 («t* + (Wg*o«o - Wb*oWb«o) + ft-c) 
+ y (^ f O_0 2 + /i.e.) + — - M^)(D_(D 2 + fc.c.) (4.10) 

+ y (^W- - ^) 2 4>_<1) 2 + ft.c.) + ff 3 (tf'<t>_(W 2 + + - W+0+W+0+) + /i.e.) + . . . , 
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Figure 2: One-loop graph with derivative vertices (denoted by filled squares), with = (pi + pif ■ 

where K denotes the non-relativistic field for the K + meson, W K = (M| - A) 1/2 , and the ellipsis 
stands for higher-order terms in e. Note that, by construction, the above Lagrangian contains 
all terms up to and including 0(e 4 ), allowed by the symmetries in the non-relativistic effective 
theory. 

5. K + -* 7>n decays: pion rescattering to one loop 

Here, we discuss terms of 0(a), generated by one-loop graphs of the type X.k X -dm, see 
Fig.|2]for a specific example. The technique used for the calculation of these loops is described 
in Sect. [2] for the case of non-derivative couplings. The presence of the latter does not change 
the procedure fundamentally, as we demonstrate for one specific example. Consider the diagram 
shown in Fig. [2] We restrict ourselves to the part of the amplitude proportional to the coupling 
H\ , which can be written as 

M H N ( Sl ,S2,s 3 ) = 2 Hl (C x + D x (si-s x )) f ^ W(X) ^ 



2w(l)2w(Q 3 - 1) (w(l) + w(Q 3 - 1) - Q°) ' 

(5.1) 

with (2^ = (pi + pif - Rewriting the numerator of the integrand as 

w(D -M n = (g°/2 - M n ) + I (w(l) + w(Q 3 - 1) - fig) + - w(Q 3 - 1)) , (5.2) 

we see that only the first term yields a non-vanishing contribution after the integration, since the 
second term after expansion in momenta leads to dimensionally regularized no-scale integrals, 
and the third term is antisymmetric with respect to 1 — > Q 3 - 1. We finally obtain 



MS 



%(S U S2, S3) = 2#i(y - M„yC x + DJs 3 - s x ))J + Js 3 ) , (5.3) 

where the function 7 + _(i 3 ) is given in Eq. (13.61 >. 

The generalization to other derivative couplings is obvious. The complete one-loop represen- 
tation for the decay amplitudes Mn and Mc, up-to-and-including terms of 0(a e 5 ), reads 

M^~ loop) (si,s 2 , s 3 ) = B N i(s3)J+-(s3) + B N2 (s3)Joo(s3) + {Bm(si,s 2 , s 3 )7+o(si) + («i <^ s 2 )} , 

M ( c~ lOOP \s u S 2 , S3) = B C l(S3)J ++ (S3) + {B C2 (SI,S2, S3)J+-(Sl) + B C 3(si)Joo(si) + (Si o s 2 )} , 

(5.4) 
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where 



B m (s 3 ) = 2[C X + D x (s 3 - s x ) + F x (s 3 - s^MHa + H x {-j- - M K ) 

*3 



+ H 2 



B N i(s3) = [Coo + A)o(>3 - soo) + foo(>3 - soo) 2 ]|g + Gi(p° - M„) 



+ G 2 (p°-M n f + G 3 f(l- 4 -^) 



i r Q° A 

B m (si,s 2 , s 3 ) = 2[C +0 + D +0 ( Sl - s +0 ) + F +0 ( Sl - 5 +0 ) 2 ]{Go + d[ -yf 1 + ) _ M - 



+ G 2 








+ G 3 




\ 2 

-rf) 




1 p ^fo) r 
3 M/f l 


- s 2 + 


* 1 



x{g! + G 2 [(l + ^)e ( » - 2M ff ] + G 3 [2p? - (l - 



+ 0(A 2 ), (5.5) 



and 



ficife) = [C ++ + D ++ (s 3 - + F ++ (i 3 - 5 ++ ) 2 ]{tf + ffi(^3 - M n ) 

+ // 2 (^-M,) 2 + i/ 3 f (1-^)}, 

Bc2(*i, s 2 , s 3 ) = 2[C + _ + D + _(ij - + F + _(i, - 5 + -) 2 ]{tf + H ^~ M * 

\(Q\ \ 2 Q 2 / 4M 2 ,] r/g ( » \ 2 Q 2 , 4M 2 , 

Baisi) = [C x + D»(J! - + - ?,) 2 ]{Go + Gi(p? - M„) 



, Q 2 4M 2 
+ G 2 ( P ?-M ;r ) 2 + G3-^(l 



*1 



)}■ 



(5.6) 



Here we have added terms of 0(e A ) (oc F,-, in a canonical extension of our notation) to the nn 
amplitudes, without having them formally introduced on the Lagrangian level. They contribute 
terms of order (s-s) 2 to the S-waves. Only some terms of 0(A 2 ) have been neglected in Eq. ( I5.5I ). 
Note in particular that there are no contributions of nn D-waves at this order. 
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Figure 3: Two topologically distinct non-relativistic two-loop graphs describing the final-state nn rescattering in the 
decay K — > 3fi, with Qf = tf x + g^. The positions of the leading Landau singularities of graph A are discussed in 
|Appendix F| 



6. K + -» 3n decays: pion rescattering to two loops 

6.1. The diagrams 

There are two topologically distinct two-loop graphs that describe pion-pion rescattering in 
the final state, see Fig. [3] In order to ease notation, we set Q 1 = + q^, s = Q 2 . Further, 



g U = ^ ^ q 2 = : *V . (6.1) 

2M K 4M 2 K 

In addition, throughout this chapter we consider the case of the non-derivative couplings only. In 
this case, in the rest frame of the kaon it is possible to express both diagrams shown in Fig. [3]in 
terms of a single variable s and one has 

M 2 -£° ps (s) = M% c (s) + M B NtC (s) . (6.2) 

The diagram in Fig. [3i, apart from the prefactor containing coupling constants, is given by a 
trivial product of two one-loop diagrams, which were already given in Eq. (13.6b . 

M B N>c (s) cc J ah (s)J cd (s) ■ (6.3) 

Obviously, in the non-relativistic framework Ai^, c (s) is therefore ultraviolet finite and of order 
e 2 . 

In the remaining part of this section we discuss the evaluation of the non-trivial contribution 
Ai^ c (s), which stems from Fig.[3j\. It is in particular shown that - up to a low-energy polynomial 
and an imaginary part which does not contribute at the accuracy we are working - c (s) is 
given by the function F(M a , Mb, M c , M^; s) defined in Eq. (16.23b . So, the reader not interested in 
the details of the derivation may directly proceed to Eq. (I6.23l >. 

6.2. Evaluation of the generic two-loop function 

The quantity Ai^ c {s) is proportional to the generic two-loop function 

AtC, ■ ' ' 11 



C d D l d D k 
J (2n) D i QMf 



i 2w a (-l - k) w a (-\ -k)-M K + f+ fc° 2w b (Y) w b (l) - 1° 

,64, 



2w c (k) w c (k) - 2w d (Q - k) Wd (Q -k)- Q° +k° ' 
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Performing the integration over the fourth components of the momenta, we can rewrite the above 
integral in the rest frame of the kaon as 



I J ab (M K -w c (k),-k) 

2w c (k)2w rf (Q - k) w c (k) + w d (Q - k) - Q 



where J a b(Lo, L) denotes the inner one-loop integral 



J ab (L ,L) = I — '— l - - . (6.6) 

1 - J 2w a (L - \)2w b (\) w fl (L - 1) + w b (l) -La 



In the calculation of the inner integral, we proceed as in the case of equal masses, see Eqs. (12.5b — 
( 12.91 ): we first shift the integration variable according to 



1„ ... . Ml- Ml 



1^1+-(1+5 L )L, 5 L = "-JT^- L 2 = (L°) 2 -L 2 . (6.7) 

Further, using the identity 

1 1 1 



v a w b I L - w a - 



4w a w b I Lq - w a -w b L° + w a + Wft LP + w a - Wb L° - w a + Wb ' 
= 2L> l 2 - (1L/L ) 2 - k\ ' (6 ' 8) 

where 



"<.= ^ 2 + (i(l-<J L )L-l)\ w,= ^M 2 + (i(l+<J L )L + l)\ ^ 



/l(L 2 ,M 2 ,M 2 ) 



4L 2 

(6.9) 

applying threshold expansion and rescaling the component l\ (see the discussion after Eq. ( 12.8b ). 
we find ^ 

J ab (L°,L) = -4= f 7t4 tAt • ( 6 - 10 ) 
2 VI 2 J ( 27r ) I 2 - k 2 Q 

Finally, performing the 1 integration, we reproduce Eq. (13.6) in the limit d — > 3. 

In the two-loop diagram Eq. ( 16.5b the components of the momentum W, corresponding to 
the one-loop subdiagram, are equal to L° = Mk - w f (k) and L = -k. One may further ensure 
that Icq = Aic(A 2 - k 2 ), where the quantities A 2 and A* are given by 



A(M 2 K ,M 2 ,(M a + M b ) 2 ) a = Ml / (M a -M b ) 2 \ 

4M 2 ' * 2(M 2 + Ml) - (M a + M b ) 2 - s k \ St I ' 



a 2 = : *' — / — — . a, = 

f2 i Jjf2 



ft = M£ + M f z - 2M K w c (k) . (6.11) 
Rescaling 1 — > (Aj) 2 I> we finally arrive at 

Jab(MK - WM, - k )=^J^ FTk 2 ^ 2 " (6 ' 12) 

Now we have to insert the above expression into Eq. (16.5) and calculate the two-loop inte- 
gral. To this end, we rewrite the denominators in the outer loop, using again Eq. ( 16.81 ) with 
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w a ,Wb — > Wd,w c and LP — > Of. The integral over the last three terms does not vanish any 
more in dimensional regularization, since these are multiplied by the inner loop, which is not 
a low-energy polynomial. It can however be easily checked that, expanding everything but the 
non-polynomial factor (l 2 + k 2 - A 2 ) -1 coming from the inner loop, one obtains a result of the 
form P(s) = P(s)(-A 2 - ie) d ~ z , where P(s) is a low-energy polynomial in the variable s with real 
coefficients and a simple pole in d - 3. Consequently, 

1 r d d l d d k Af 2 ~' 1 1 _ 

M(S) ~2^J W l 2 + k 2 - A 2 p2 - (pQ/gO) 2 - q\ + P(S) = M(S) + P(S) ' 

(6.13) 

where k = p + ± (1 + 5)Q, and 

, A{s,M 2 c ,M 2 ,) M 2 -M 2 
q 2 o= i-^' -• (6-14) 

6.3. Renormalization 

At this stage, it is appropriate to discuss the freedom in the definition of M{s). For example, 
one may add an arbitrary low-energy polynomial of s with real coefficients to M{s) - this would 
amount to a renormalization of the K — > 3n vertices G,, Hi in the non-relativistic effective 



Lagrangian. One may use this freedom to remove the real part of the polynomial P(s), which is 
of order a 2 and which contains an ultraviolet pole at d — 30 On the contrary, the imaginary part 
can not be removed in this manner. So if the imaginary part of P(s) were divergent at d — > 3, 
it would constitute a major problem for the validity of the framework. However, this does not 
happen: as can be easily seen, the imaginary part is ultraviolet-finite. Moreover, at the accuracy 
we are working, one may neglect this imaginary part altogether, because its contribution to the 
decay rate starts at order a 3 , that is beyond the scope of the present paper. 

To summarize, the two-loop function Ai(s) at the accuracy we are working can be replaced 
everywhere by 

F(M a ,Mb,M c ,M d ;s) = M(s)-ReM(s t ), (6.15) 

where s, = (M c + Mj) 2 . Here, we used the above-mentioned freedom to normalize the real 
part of F(M„, Mt, M c , Mj; s) to zero at s — s,. The difference between F and M is a low- 
energy polynomial. The real part of this polynomial can be removed by renormalization, and the 
imaginary part is ultraviolet-finite and does not contribute at the required precision. 

6.4. Integral representation 

The rest of this section is devoted to the calculation of the function F. After shifting the 
integration variable in Eq. (16.13b according to 

we arrive at 

M) = j_ r_d^i_d^k_ i n(x) 

2 VI J (2nY (2nY (k 2 - ? 2 ) (P + 0^! Q 2 + , _ A 2) ' 



4 Renormalization of the couplings Gj, Hj first occurs at two-loop order in our framework. 
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where 



x = k2 + (kQ^ + kQ e „ (1 + ^ 

S yfs 



N(x) 



(Ski*)) 



2 y/s k (x) 



9 9 / 9 (l+<5) 2 9 \' /2 

^(pc) = M K + Ml - 2M K [M 2 C + Q 2 + xj . 



(6.18) 



In order to evaluate the integral in Eq. (16.17K we expand the numerator N(x) in the variable x 
and integrate term by term. Using Feynman parameterization to combine the two denominators, 
we obtain 



1 ^ 1 d" 
M(s) = — — > — — N(x) 
2 V? U nl dx " 



Jn{s) : 



r=0 



Us) = I dyy ~ d 'i l + T")~ ^ ->Ctfy. *) - 

where g(y, s) is defined by 



g(y, S ) = -(l-y>7 2 -yA 2 + ± 



h(i-y)Q 2 (i+s) 2 



1 + 



(6.19) 



(6.20) 



and the first few coefficients in the expansion are given by 
f d (y,s) = l, f d ( y, s) = w -^(l + —) + y, 



(6.21) 



(2 + d)g z (y, s)(, f 2Q z a _ 3QV 
4(l-fl?) 



1 + 



jV \ dg(y,s) ( 
ds ' d(d + 2)s 2 ) 2{\-d)\ 1 



Q 2 (2ay+/3 2 s)\ 



Here, 



1-y 



1 + 



Q 2 (G ) 2 „ . o2 i -K 1 -?) 

(i + *?) a 



2s 



■y(l +6) 



. (6.22) 



In order to eliminate the singularity at y = when d — > 3, we first integrate by parts and omit 
the surface term, which is a low-energy polynomial in s. One may further verify that in the 
limit d — > 3, to a low-energy polynomial, the function F is given by the following integral 
representation 

F(M a ,M b ,M c ,M d ;s)= \ f Hy, s)(ln g(y, s) - In g(y, s,)) . (6.23) 
256n 3 yfs Jo y/y v y 

The function !F(y, s) is given by an infinite sum 



T(y, s) = ) Tniy, s) , T„(y, s) = — —N(x) — = 0( 

«! dx- ,=o dy\ > 



e 2 " +2 ) , (6.24) 
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Figure 4: Successive approximations for Re F(M n , M„, M n , M n ; s) = ReF(.s) in the equal-mass case. We show the 
leading-order (maximum n = 0, dashed), next-to-leading-order (n = 1, dash-dotted), and next-to-next-to-leading-order 
(n = 2, full line) approximations, see Eq. 16.24) . 



and f^ n \y, s) = fj (y, s )\ d -y I n particular, for n — we find 



U=3 

^ 2 (s ,M 2 a ,M 2 b ) 1 d{ g(y,s) \ 

To(y, s) = z — , (6.25) 

*° _ ii^li Q 2 ^i + zQ: 7 



where 



5 = M| +M c 2 -IMk^MI + ^(1 +5) 2 . (6.26) 

The series in n for the function F(M a , Mb, M c , M^, s) converges rapidly: we display the succes- 
sive approximations with maximum n — 0, 1,2 for the equal-mass case in Fig. |4] It is seen that 
there is almost no difference between the two approximations with maximum n — 1 and n — 2. 

The function F(M a , M b , M c , M d ; s) given by Eq. (16.231 1 starts at 0(e 2 ) and thus does not 
violate power counting. Note also that, expanding the function T in powers of — = 0(e 2 ), it is 
possible to perform the integrals in y analytically at each order. The result at next-to-leading order 



is given in Appendix B Analytic properties of this function are considered in Appendix C and 



the comparison to the relativistic approach is discussed in Appendix G 



6.5. Threshold behavior 

Finally, we display the singularity structure of the function F{M a , M/,, M c , Mj; s) near the 
cusp. It can be shown that 



Lt 1/2 (s, M 2 , M 2 ) iA 1/2 (s 0t , M 2 Mh 
F(M a , M b , M e , M d ; s) = — d - — -2 b - + O(q 2 ) , (6.27) 

lOTTS lOTTSQt 



where so, denotes the function s^(s) in Eq. (16.26b . evaluated at s — s,. Hence in the vicinity of the 
cusp, the two-loop diagram is given as a product of two factors: the first factor describes the cusp 
emerging in the outer loop, while the second factor is the inner loop evaluated at the threshold 
s = 4M 2 . Thus, the above two-loop diagram satisfies the threshold theorem ill l2lll22ll23l 12411 . 
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Figure 5: A non-relativistic two-loop graph describing the decay K + — » n°n°n + that involves a six-particle vertex. The 
intermediate state contains the particles n°(k[ ), 7T°(k2), and n + (Pn —k\— kz). 



7. Three-pion couplings 

Pion rescattering in the final state also contains a contribution from diagrams of the type 
shown in Fig. [5] These diagrams contain a vertex describing the interaction of six pions at the 
same point. In the non-relativistic theory, such interactions are generated by a Lagrangian of the 
type 

iF ot(Oj) 2 O + <I>2 + i. 

where the ellipsis stands for terms with space derivatives. 

For demonstration, let us consider the diagram with a 7T°n°7T + intermediate state, depicted in 
Fig. |5] The contribution of this diagram to the decay amplitude Mn is a constant, which in the 
center-of-mass frame of the kaon Pk = (Mr, 0) is given by 



£3* = jF Ol(<Dj) 2 (D + (l)2 + 7 ^<l)l((Dl) 2 cD_(D 2 + + ... , (7.1) 



war, v 1 r F f d% f d d k 2 1 



1 



i) d 2wo(ki) 2w (k 2 ) 

(7.2, 



2w(ki + k 2 ) wo(ki) + wo(k 2 ) + w(ki + k 2 ) - M K 

Obviously, the diagram is a constant. Its real part, which is in fact divergent for d — > 3, can 
therefore be absorbed in a redefinition of the coupling constant Go, and will not be considered 
further. We concentrate on the finite imaginary part, for which we may take d — > 3 and rewrite 
the above expression in the relativistically covariant manner 

. ^ 3 

Im M%(si , s 2 , s 2 ) = G F J f] d% e(k?) 6(kf - M]) 6 (4 \P k -k x -k 2 - h) 



= 1 

2 

+ 0{K,{M K -3M n y). (7.3) 



768 V3tt 2 



Similar to Ref . [21], for the rough estimate of the magnitude of the coupling Fq we use matching 
to chiral perturbation theory. At lowest order this results in 

Ml 

F = + . . . , (7.4) 

where the ellipsis stands for higher chiral orders, and F„ = 92.2 MeV is the pion decay constant. 
Since also the imaginary part of the diagram in Fig.|5]is a constant, it is possible to eliminate this 
contribution by allowing the constant Go to have a small imaginary part 

ImGo (Mk-3M,) 2 M 2 t c 5 



Re Go 768 V3V 8F* 



^1.5-10 . (7.5) 
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Note that Eq. ( 17.5b does not contain effects from the n + n + n~ intermediate state, which would 
contribute a similar term oc Hq /Go- As the effect turns out to be so small, we will neglect it below 
and continue to assume that the couplings in the non-relativistic kaon effective Lagrangian are 
real. Note that our result is in qualitative agreement with Eq. (80) ofRef. [21], although the exact 
coefficients differ. 

We wish to stress that the non-relativistic theory without six-particle couplings is in general 
not self-consistent (see, e.g., Ref. jjj}] ). The reason for this is that there exist divergent two-loop 
graphs for the process 3n — » 3n that require the introduction of a six-particle counterterm in the 
Lagrangian. These diagrams arise, however, at order a 4 in the expansion in (small) nn scattering 
lengths and thus are beyond the scope of the present article. 



8. K + -» 3n decay amplitudes to two loops: final result 

In order to find the representation of the decay amplitudes at 0(a 2 ), we draw all possible 
graphs for K + — > ji n Q n + and K + — > 7r + 7r + 7r~ at one- and two-loop order, comprising the topolo- 
gies discussed in detail in Sects. |5]and|6l for an explicit display of the different pion intermediate 



states, see Refs. l22l 12311 . Further, in the graphs of type Fig.[3]\ we retain only the non-analytic 
piece F(M a , Mb, M c , Mj; s), whereas the (infinite) polynomial is included in the tree-level cou- 
plings G,, Hj. The choice of the particular representation Eq. (16.231 is equivalent to the choice 
of a renormalization prescription. Note also that the resulting modification of G,, //, is of order 
a 2 , so that if one uses the modified couplings also to calculate one- and two-loop diagrams, the 
results change at 0(a 3 ) and 0(a A ), respectively, that is beyond the accuracy we consider. 
The amplitudes up to and including two loops are given by 

M 00+ - MT + <~ lo ° P + M 2 - loops + ... [K + ^ ttW] , 

M ++ - = M£ ee + M l c loop + M 2 c loops + ... [K + -» ttVtT] . (8.1) 

Our amplitudes are normalized such that the decay rates are given by 

*'^ 2 " mp <- p > me ni0Fr (82) 

The tree-level contribution to the K + — > n°n°n + and K + — > n + n + n" decay amplitudes up to 
(9(e 4 ) is given by Eq. ( 14.81 1. The one-loop amplitude fully consistent up to 0{ae s ) is displayed 
in Eqs. ( 15.4b . ( 15.5b . and (15.6b . Finally, the two-loop contributions to the amplitude, which are 



by far the lengthiest and most complicated part, are given explicitly in Appendix A These 



are complete up to 0(a 2 e 4 ), but contain partial terms of 0(a 2 e b ) and 0{a 2 e*), retained for the 



following reason. As we have seen in Sect. 13.21 no reasonable description of the nn S-wave 
scattering amplitudes is possible without including the effective ranges. We have therefore kept 
all linear energy dependences in the S-wave interactions, generating in combination terms up to 
0(a 2 e & ). In the two-loop graphs, we neglect however terms induced by couplings quadratic in 
energy (oc G2/3, H2/3, nn shape parameters) as well as higher orders in P- (or even D-)waves. 

Prior to using this representation in the fit, one eliminates the coupling constants C,, D,-, in 
favor of the nn threshold parameters through the matching conditions discussed in Sect. [3] This 
representation depends on the eight real K — > 3n coupling constants H„ G, and the threshold pa- 
rameters for nn scattering, which in the end ought to be determined from a fit to the experimental 
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data. It is this amplitude, fully documented here for the first time, that, amended with the radia- 
tive corrections provided in Ref. 112411 . has been used in the analysis of the NA48/2 K + — > n°n°n + 
data lfl7ll . which in combination with results from K e 4 decays led to the determination of nn 
scattering lengths quoted in Eq. ( 11.21 . 



9. Comparison with other approaches 



The low-energy expansion proposed here is closely related to early work performed in the 
1960s by many authors, who used S -matrix methods to investigate the production of particles 
in the threshold region, in particular also in K — > 3n decays. We refer the interested reader to 
the article by Anisovich and Ansel' m l45ll for a review. The framework presented here may be 
considered an effective field theory realization of these approaches. The most notable differences 
are: 

1 . In these early approaches, the occurrence of a cusp in the decay rates went unnoticed. 

2. It is explicitly stated in Ref. [45] that there are no leading Landau singularities in the graphs 
of the type Fig.[3j\, which does not agree with what we find. As a result, our expressions 
for the two-loop integrals do not agree with the ones presented in these early articles. 

3. Our approach allows for a straightforward inclusion of the effects of real and virtual pho- 
tons, see Ref. i24ll for the actual evaluation of the pertinent matrix elements, and Ref. lfl7ll 
for applications in the data analysis of K — > 3n decays. 

A comparison of the present framework with the seminal articles of Cabibbo and Isidori J2, 
I2H1 was already provided in Ref. [22], and we refer to that article for details. Here we only recall 
that the decomposition of the decay amplitudes J2,|2l|] into the form M = Mq + v C d(s)Mi, with 
Mo,i holomorphic, is valid near the threshold only, as is detailed in Appendix H Furthermore, 



the leading Landau singularities found here were missed in Ref. 12j,|21|]. However, the analysis 
in Ref. [17] shows that, at the accuracy considered, these effects are immaterial for the main 
problem at hand, the extraction of the nn scattering lengths from the cusp strength. Various 
aspects of K — > 3n decays are presented in the recent articles Refs. H Ey4i|5a|5il|5l|5l . 
A comparison of these works with the present approach is provided in Refs. H22L I2.3L 12411 ._see also 



Ref. Ill for a careful comparison of the results of Refs. I150l Dll 15211 with those of Ref. 12411 . We 
refer the interested reader to these articles for details. 

A dispersive approach has recently been undertaken to K — > 3n (and 77 — > 3n) decays 
and an analytical expression for the decay amplitudes for Ki — > 3n valid up to two loops in chiral 
perturbation theory is already available. On the other hand, it seems to us that the effects of real 
and virtual photons are very difficult to incorporate without a Lagrangian framework. The main 
difference to the framework presented here is the fact that quark mass effects are treated in a 
perturbative manner in Refs. 115 411 . 



10. Summary and conclusions 

We summarize the main findings of our investigations as follows. 

1. nn rescattering generates a cusp in the decay distribution of K* 1 — > n ^^ decays and 
allows one to measure nn S-wave scattering lengths «o and A]. This cusp has been 
investigated in recent years by the NA48/2 collaboration in great detail JHH^, 171. and the 
two scattering lengths were determined in this manner. 
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The analysis of the NA48/2 collaboration relies on two theoretical frameworks 
which relate the decay distribution to the scattering lengths. In this article, we spell out 
details of the method that underlies the matrix elements worked out in Refs. I22I l23ll and 
used in Ref. [ 17] for the data analysis. 

We investigate the decays K — > 3n within an approach which is based on non-relativistic 
effective Lagrangians. It enables one to systematically calculate the decay amplitudes 
K + — > 7r () 7r°7r + and K + — > n + n + n~ in a double expansion in small momenta of the decay 
products and in the threshold parameters of the jm scattering (scattering lengths, etc.). Fur- 



thermore, it allows one to take into account radiative corrections in a standard manner [24]. 

4. Because it is a Lagrangian framework, the strictures of unitarity, analyticity and of the 
cluster decomposition are automatically taken into account. Furthermore, the amplitudes 
contain the Landau singularities automatically. 

5. In this article, we provide the amplitudes for /T ± — > 7r°7r 7r ± and for K* — > t&t^ti' in 
closed form, at two-loop order, up to and including terms of order a 2 e 4 , with the most 
important terms of order a 2 e 6 and a 2 e & also retained. The extension to Kl — * 3n decays is 
straightforward. 

In conclusion, the framework presented here is very suitable to investigate the beautiful set of 
data on K — > 3n decays collected by the NA48/2 collaboration - it has all the bells and whistles 
provided by quantum field theory. We expect that the empirical parameterization of the Dalitz 
plot distribution published recently [56] would even allow one to study the convergence of the 
a, e expansion and to narrow down the emerging theoretical uncertainty in the scattering lengths 
«o,2- This endeavour is, however, beyond the scope of this work. 
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Appendix A. K + -> 3n amplitudes up to 0(a e ) 

We have calculated the two-loop contributions to the decay amplitudes for K + — > 7r°7r°7r + and 
K + — > t: + 7: + 7t~ including derivative couplings, along the lines explained in detail in Sect. [6] The 
results we find, complete at order a 2 e 4 and containing additional terms of 0(a 2 e 6 ), 0(a 2 e & ) (see 
the comments in Sect. |8), read 

M^ ops = M^ c (s u s 2 , s 3 ) + M^cisu s 2 , s 3 ) , (A.l) 

where 

M A N = 4H'_;'_(si)C + -(5 + 3 -)CAs3)F + (M„, M„, M„, M n ; s 3 ) 

1 Q 2 

- 4[-// 1 C + _(s 3 + -) + 2M K H';'_(s 3 )D + .]c x (s 3 )-^ F + ( M «> M *> M *< M - ft) 

Q 4 

+ AMkH^C^)— |- ff >(M„, M„, M„, M„; s 3 ) 
(t; 3 ) 

+ 2G"( i3 )C,(S 3 + -)C,( i 3)^ + (M ;r o, M,o, M„ M*; * 3 ) 

O 2 

+ 2[GiC,(St) - 2MffG"(s 3 )AjC JC (j 3 )-rj F^{M^,M^,M n , M„; s 3 ) 

Q 3 



- 4M Jf G I D x C x (s 3 )-^- Mj, M„ M„; s 3 ) 



2° 3 ) 2 + 

+ 2//"( 53 )C ++ (S 3 f -)C x ( i3 ) J F + (M„ M„, M„, M„; s 3 ) 



Q 2 

+ 2[H l C ++ (s + } -) - 2M K H"(s 3 )D ++ ]C x (s 3 )-^ F ( ^\m k , M n , M n , M n \ s 3 ) 

y 3 



- AM K H x D + +C x {s 3 )—2- F {2) (M n , M„, M„, M n ; s 3 ) 

\H 3 ) 

+ 4G;',' 1 (.9 3 )C + o(S ( 3 lo )Coo(5 3 )F + (M ?r , M^o, M#, Mj>; s 3 ) 

1 Q 2 

- 4[-GiC +0 (^°) + 2M K G%(s 3 )D +0 ]Coo(s 3 )-l F^(M n , M^,M^,Mj>; s 3 ) 



Q 



+ 4M*G 1j D + oCoo(s 3 )-^ Ff\M„, M„ , M n o , Mj> ; s 3 ) 

+ {4i/;'!( i |)C,(J[ )[C +0 (ii) + «2, s 3 )]F (M n , M„, M n , M^sx) 

dHiC x (Sf) + 2M K H':'_(4)Dx)c +Q (si& + 2i/;'!(4)C ;c (s; )£ + o(si, s 2 , ft) 
2 Si 

o 4 

+ AM K H\ D x C+q(s i ) ^ 2) (M„ M„, M ff , M,o; j,) 
+ 4G^(^)C + o(S? + )[C + o(ii)-£;o(^i^2^3)] J Fo(M,o,M„M^,M,;i 1 ) 
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r 1 O 2 
- 4 (-G!C +0 (S? + ) + 2M K G'^{ S - x )D + ^C +0 ( Sl )-± - 2G^(sr)C +0 (s° + )E +0 (s 1 , s 2 , s 3 ) 



o 4 

+ 4M*G 1 D +0 C +0 (s 1 )— ^- FfiM^Mn, M*,M*; si) 
+ 2G"(4)C o(s| )[C + o(si) + Ei,(si,ii, ss^M^.M^.M^M^; si) 



+ 2 



o 2 

(GiCooOf) - 2M jC G"«)Doo)C +0 (si)^ - 2G"(s+)Coo(s|°)£ + o(si, s 2 , s 3 ) 
xF^M^.M^M^M^si) 
4M J ,GiD 00 C + o(si)-^ F { q\M k o , Mjfi , M„, ; s x ) + (sj ~ s 2 )} , (A.2) 



At* = 4//'(s3)C,(s 3 )C + _(s 3 )7 2 _(^) + G(s 3 )C 00 (s 3 ) 2 ./ 2 (s 3 ) + 2G(s 3 )C,(s 3 ) 2 y o(^)A-(s3) 
+ 2H'(s 3 )C x (s3)C o(s 3 )J + -(s 3 )J o(s 3 ) + {4G'(s 1 )C +0 (si) 2 7 2 (s 1 ) + (sj «-> s 2 )} , (A3) 



Al£ = 2G"( S3 )C x (t 3 + )C ++ (s 3 )F4M jfi ,M jfi ,M n ,M K ; s 3 ) 



+ 2[G!C x (sJ + ) - 2M*G"(s 3 )Z)JC ++ (s3)-f F^M^o, Mj>,M n , M n ; s 3 ) 



(G3) 2 



■ 4M*G 1 A c C ++ (s 3 )— ^ F & \M n o,M„o,M„ M n ; s 3 ) 



+ 4#;':(s 3 )C + _(s 3 ++ )C ++ (s 3 )F_(M, r , M„, Af ff) M„; s 3 ) 

1 Q 2 

- 4[-// 1 C + _(S 3 ++ ) + 2Mk//;"(s 3 )D + _]c ++ (s 3 )-| F ( l\M n ,M n ,M n ,M n - s 3 ) 
z tJ 3 



+ 4M jC // 1 D + _C ++ (s 3 )— ^- F_ (M„, M„, M„, M n ; s 3 ) 
(G 3 ) 

+ {4//; , :(s 1 )C + _(s|-)[C + _(s 1 ) - F + _(s 1? s 2 , s^F+^M^M^M^ s x ) 
-4 



1 o 

(-//!C + _(s|-) + 2M J ,//;'l(s 1 )D + _)c + _(s 1 )^ 

"1 



■2//;':(s 1 )C + _(s[-)£ + _(si,s 2 ,s 3 ) 



^ 1) (Af„Af jr ,Af„Af,;j 1 ) 



+ 4M*# 1 D + _C + _(s 1 )— ^- Ft\M n , M„, M K , M n \ s x ) 
(Gi) 

+ 2G"(si)C^(s|-)[C + _(s!) - F + _(s x , s 2) s 3 )]F + (M^ , M„o , M„, M n ; s x ) 



+ 2 



(GiC(ST) - 2M jr G"(s 1 )D,)C + _(s 1 )-i 
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+ 2G"(s l )C x (s + l -)E + -(s u s 2 ,S3) 



(G?) 2 



4M K GiD x C + -( Sl )-^- Ff\M„o,M n o,M„ M„; si) 



+ 2i/"( il )C ++ (^ + )[C + -(^i) + £+-(si, s 2 , s 3 )]F + (M„, M„, M„, M n ; s x ) 



+ 2 



(i/!C ++ (^ + ) - 2M K H"( Sl )D ++ )c + -( Sl )j± 

Gi 



2H"( Sl )C+ + {s- x + )E + -(s x ,S2, s 3 )Y ( l\M n , M„, M„, M„; si) 



- 4M K HiD ++ C + -(Si)-^- F+\M n , M„, M n , M n \ Sl ) 
+ 4G^(si)C + o(S?°)C,(si)F + (M„,M a o,M a o,M B o; 5i) 



1 o 2 

- 4[-G 1 C +0 (S? ) + 2M*G + ' (5i)Z)+o]c ;c (5i)--£ F+^Af*, , A^ , M n * ; Sl ) 
z Gi 



(G?) 2 



+ 4M jr GiD +0 C J ( Sl )-^- Ff{M K ,Mj,Mf,M#\ Si ) + (si «-» s 2 ) 



(A.4) 



Alp = //(^)C ++ ( i3 ) 2 4 + (^) + {4//'(5i)C + _(50 2 ./+_(5i) + G( il )C,( il )C o(^i)4)(^) 
+ l[H'{si)C x {si) + G(si)C + _(si)]c I (ji)7 + _(ji)7oo(Ji) + (*i <-» S2)) ■ 
We have used the following abbreviations: 



M K 



C n (Si) = C n + D n ( Si - s n ) , sf = M z c + Mf - Si + -±( Si + 2Qf - M z c + M z d ) , 

£+-(51,52, 53) = £+- Sl ^H 4p , £+0(51,52, 53) = £+0 



£+ Cs , l,52,5 3 ) = £+0 



(51 ± A n )(s 3 - 5 2 - AJ 



5l(53 - 5 2 - A„) 



2M K Q\ 



+ A, 



2M^G 



+ A^ 



2 O 2 4Af 2 
G(5,) = G + Gi (p° - M„) + G 2 (p? - M„) + G 3 -f(l , 



AMI 



(A.5) 



fffa) = H + Hi (p° - Af„) + ff 2 - M„) + #3^-(l -) , 

G'( Si ) = Go + Gi(^(l + ^) - M n ) + G 2 [(^(l + ^) - M,) 2 + ^ X(s„M-. M;. ) 



+ G 3 



G° x r/G° x2 Q 2 
f-M n ) + H 2 [(^--M^ ■ ' 

(l-p?f4(.-^)l 



ff'(5,) = tf + - M„) + H 2 [(f- - M,) 2 + ^. (l - _JL) 
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G"( Si ) = Go + Gi(^ - M„) , G£fo) = G + G^-(M K - M a - M b ) - Jj-} 



4fi 



G"(st) = Go + Gi(^^ - M„) , GZ(sj) = G„ + - M a - M 6 ) - , 



H"{ Si ) =H + - M„) , = Ho + H^-(M K - M a - M h ) - ^L) , 



4gy 



KM) = H o + H^(M K - M a - Mb) - S -j^) . 



(A.6) 



We remark that the representations of the "double bubbles" At^ c are even complete to (9(a 2 e 6 ) 
if the polynomials C„(i,) are amended by shape parameter terms . . . + F„(sj - s n ) 2 : P-wave 
contributions only start at 0(a 2 e^). In contrast, in the "genuine" two-loop graphs At^/ C , terms 
oc E n F < - 2) have not been included, although they contribute in principle at 0(a 2 e 6 ). We remark that 
the two-loop amplitudes for the Kl — > 3n decay channels at the same accuracy can be retrieved 
from Ref. [29], where the 77 — > 37r amplitudes are discussed within the same formalism, with the 
obvious replacements as described in Ref. ll23ll . 

Fi(. . . ; s), F?\. ..;s), F?\. ..;s) stand for the integrals F(. . . ; s), F m (. ..;s), F {2) {. ..;s), 
evaluated at Q 2 = AiM^, M 2 , s)/4M 2 l , with i = ±, 0. The first few terms in an expansion in e of 
the analytic expression for these two-loop functions are displayed in Appendix B 



Appendix B. The loop functions F, F (1) , and F {2) 



The functions Fj(. . . ; s), F^(, . . ; s), . . ; s) introduced in |Appendix A| are explicitly 
given by 



7 (2), 



F(M a , M b , M c , M d , s) = N 



3Q 2 2 
2v/i +pfo- JfaiPfi ~ 2 1ofo) 



+ K(X 3 f 3 + X 2 f 2 + X\f\ + X f ) 



+ 0(e b ) 



F ( - 2 \M a ,M b ,M c ,Md,s) = 'j 



N_ 

10 
N 



F m (M a , M b , M c , M d , s) = + <5)[(10v - p)/j + (5p + 2q 2 )f ] + 0(e 4 ) . 



-^(2v 2 / 3 + 3vpf 2 + {p 2 - 2vq 2 Q )f l - pqlfo) 
(v/ 3 + (p- 2v)h + (4v - 2p - + 2(p + q 2 )f ) 



(B.l) 



with 



7C = 



1 ^ 2 Oo,M 2 ,M?) 



2567r^ 5o ^ A2 _(i^ Q2 ' 



2(M 2 + M 2 ) - (M a + Mb) 2 -s s - (M a - M,,) 2 s J s -M\- M 2 



Mi 



fa = 4(vi + v 2 - v 2 + h) , /1 = - (yi(vi - 1) + y 2 (v2 - 1) - - 1) + ft) , 
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■ 1 (3p/l - ^/ ) , /3 = ^ [3(V^ + p 2 )/! - pg2/ ] 
1 / 1 + Q 2 /* \ _ 2 _ ^l- ?3> S ) A 2 



fe=rN t y. , Q 2 = T^T^' Q = Q > s, = (M c + M d > 

2 \l+Q 2 /i,/ 4M 2 

, — 1 , _ /— — 1 -p + 

Vi = y-yi arctan , z = 1,2; V2 = y— ?2 arctan _ , yi^ = — 



y< V-3'2 



v = -— (M 2 + A 2 ), p = <7 2 -A 2 + ^-M 2 , D=p 2 +4v ? 2 , 



X = (R-Hp)q 2 , X l =H(p 2 -2vq 2 )-R(2p + q 2 ), X 2 = 3Hvp - 3R(v - |) , 

3 / q2^ O 2 2 
X 3 =2v(Hv + R), ff = --(l + |-j, R = ^f (1+S)2 > 



Q 2 (l+5) 2 \ 1/2 2 Ms,M 2 c ,M]) 
% = 



4s 



s = M 2 + M 2 - 2M k (m 2 c + 4 + ) 

, A(M 2 M 2 (M a + M b ) 2 ) Ml - M 2 , 
A 2 = — ^ — — , 6 = — , (B.2) 

4M 2 s 

compare also Fig.[3]\. The arctan is understood to be evaluated according to 

1 , 1 + ix 

arctan x = — In , (B.3) 

2i 1 — ix 

and s is given a small positive imaginary part in all arguments, s — > s + ie. 



Appendix C. Holomorphic properties of F 



The loop functions F, F {1 \ and F (2) in |Appendix A are generated by the two-loop graphs 
displayed in Fig. [3]\. The explicit expressions in Appendix B are valid on the upper rim of 
the real s-axis. Here, we show how to analytically continue them to the whole non-relativistic 
region. The holomorphic properties of these loop functions play a role in connection with the 
decomposition of the amplitude as proposed by Cabibbo and Isidori 12111 . see also the discussion 
in |Appendix H| 

The region of holomorphicity could be established directly from the explicit lowest-order 
expressions given in Appendix B - higher-order terms in e emerge from the Taylor expansion 



of low-energy polynomials and do not change the analytic structure. Here, we follow a different 
path and derive the singularity structure from the integral representation Eq. (16.231 . To ease 
notation, we set in this section 



F(s) = F(M a ,M b ,M c ,M d ;s). 



(CI) 



Appendix C.l. Holomorphic properties from the integral representation 

First we note that the prefactor T(y, s) and the function g(y, s) in the integral representation 
Eq. (16.231 1 are low-energy polynomials in [s - (M c + M,/) 2 ) and in y. Because the argument of the 
logarithm does not vanish for < y < 1 when s approaches the real axis from above at Re s > s,, 
F(s) is analytic in the part of the non-relativistic region located in the upper half of the complex 
i-plane. To analyze the singularities that may occur during the continuation to the lower rim of 
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the real axis, we note that g(y, s) has two zeros in the low-energy region, which we denote by 
yi(s) and y 2 (s), 

8(y, s) = (y - yi(*))Cy - )*(*)) . (C2) 
1 + 

where yi^C*) and v are defined in Eq. (IB. 2b [the sign convention is chosen so that yi(s) = at 
5 = 5, = (M c + M d ) 2 ]. During the analytic continuation, singularities may occur at the following 
values s of s. 

i) yi,2(S) = or 1. This may generate endpoint singularities. 

ii) If the integration contour runs between yi(s) and yiis), a pinch singularity may occur 
whenever y\(s) = y2(s). 

iii) If | lim^jyiCi)! — > oo or | lim 5 ^jy2(i)| — > °°, dragging the contour along, the integral in 
Eq. (16.23b may diverge at s and lead to a singularity in F(s). 

We now discuss these possibilities in turn. 

Appendix C.l.l. Endpoint singularity 

The quantity y\(s) [yiis)] vanishes linearly [tends to a constant] as s — > s,, as a result of 
which the derivative dF(s)/ds diverges at threshold. We conclude that F (s) is singular at s = s,. 
The case y,(5) = 1 does not occur in the low-energy region. 

Appendix C.1.2. Pinch singularity 

Concerning point ii), we note that D(s) = is a necessary and sufficient condition for the 
equality y\(s) = y2(s) to hold, see Eq. ( IB. 21 ). The quantity D(s) may be factorized in the following 
manner, 

_ P(s, q 2 v M 2 K - M at M b , M c , M d )P L (s, q 2 3 , M 2 ,; M a , M„, M c , M d ) 

6AM 2 K M 2 d s 2 

Here, P and Pi denote two polynomials of second order in s. Whereas P has no roots in the 
low-energy region, the roots of Pi, explicitly given by 

P L = A(s, M 2 , M 2 c )A{q\, {M a + M h ) 2 , M 2 ) - Q 2 (s, q\, M\; M a , M b , M c , M d ) , 
Q = 2M 2 (M 2 K -s-q 2 ) + (q 2 - (M a + M b f + M 2 )(s + M 2 - M 2 ) , (C.4) 

denote the positions of the potential leading Landau singularities of the relativistic two-loop 
graph, obtained from Fig. |3K by replacing all propagators with the relativistic ones ll57ll . We 
refer the interested reader to |Appendix F for a discussion of this point. We conclude that the 



positions of the potential pinch singularities in the non-relativistic graph Fig. |3]\ are identical 
to the location of the leading Landau singularities of the corresponding relativistic Feynman 
diagram. 

The relevant roots of the equation D(s) = can be classified according to the values of the 
masses M a , Mb, M c , M d , and M3 = (q 2 ) 1 ^ 2 . To simplify the discussion, we consider the following 
three cases, visualized in Fig.|6}\-C, see also Fig. [3] 

1. Fig. UK: 

The equal-mass case with M a - Mb = M c - M d = M3 = M„. Here, D(s) has a second- 
order zero on the real axis, at s = s\ = (M 2 , - M 2 )/2. 
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Figure 6: Two-loop graphs contributing to K + — > n + n°n° . Their holomorphic properties are discussed in |Appendix C| 
^Appendix D| and |Appendix E[ in the case where the propagators are the non-relativistic (the relativistic) ones. 

Im s 

^ + * , Re, 

v s t *s-te s p 

Figure 7: The path in the complex s-plane from the upper rim of the cut to the lower rim. The physical function at real 
values of s coincides with the boundary value of the analytic function F(s) on the upper rim. 



2. Fig. EE: 

M a - M b = M c - M„, M d = M-i- M,fi. Here, D(s) has two real roots si_ and s\_+, 

s l± = \ [Ml + 2Ml - 3Ml] ± [(Ml - M^AiMlMlAMl)]" 2 . (C.5) 

3. Fig. [S:: 

M a = M b = Mjfi, M c = M d = M 3 = M„. Here, D(s) has complex-conjugated roots at 
s — s a , s* in the complex s-plane, with 

*a = ~ K + 3M - ~ AMl A ~ 2^ V M « - Ml a )A(M 2 K ,MlAMl )\ 12 . (C.6) 

Note that the roots differ from each other by terms of 0(e 2 ). 

Let us now continue F(s) to the lower rim of the real axis, along a generic path near threshold, 
displayed in Fig. [7] It represents a circle with center at s = s, and a radius \s - s,\. We assume 
that | s - s t \/(4Ml) is small, such that the zeros of D(s) are avoided. The variable s travels along 
this circle from the upper rim of the positive real s-axis to its lower rim. 

The roots y\(s) and y2(s) then move in the complex plane along the trajectories shown in 
Fig- El When arg(s - s,) becomes equal to n, the rootyi(s) hits the original integration contour 
and enforces its deformation. When s approaches the real axis from below, the deformed path 
looks as follows (see Fig. [8): it starts at y = 0, goes in the negative direction, encircles the 
singularity at y = j ml „(s) = y\(s) and goes back to y = 1. To illustrate what happens when a zero 
of D is encountered, we consider M a = Mb = M n o, M c = M d = M3 = M„ (case [3] above), and 
envisage analytic continuation of F{s) from s* to s a along the segment of a circle with radius 
\s, - s a \. The trajectories of ^1,2(5) for the variable s moving along this segment are shown in 
Fig. [9] Since yi(s a ) = yi{s a ) (and equally for s* — > s a ), the trajectories start at the same point 
and converge to the same point - the potential generation of a pinch singularity at s = s a is 
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0.8 




-2 -1 

Rey 



Figure 8: The trajectories of y = yi(s) and y = y2(s) in the complex y-plane, when the variable s travels around the 
generic circle shown in Fig.0 For Im s < the integration contour has to be deformed as shown, in order to circumvent 
the singularity at y = y\ (s). 




Fi gure 9: The deformation of the integration contour in the variable v, as s approaches the anomalous threshold s a . 



clearly visible. The cases Q] and [2] look similar. In order to investigate whether these anomalous 
thresholds really do generate singularities in F(s), we consider the function H{s) defined by 



H(s) 



± C 1 dy_ 

Jo Vy 



dy 



g(y, s) 



In g(y, s) , g(y, s) = (y - yi)(y - y 2 ) , 



(C7) 



cf. the integral representation of F (s) in Eq. ( I6.23l l. H(s) differs from F(s) only by terms which 
are irrelevant as far as holomorphic properties are concerned, but is much simpler to discuss. We 
denote by H + {HJ) the function H(s), evaluated at the upper (lower) rim of the real s-axis, by 
continuation according to Fig. [7] We perform a partial integration and obtain 



An 

H-{s) = H + (s) + — ypf x [yi + 3y 2 ] , s > s, . 



(C8) 



Inserting the expressions Eq. (1B.21 for y^, it is seen that the second term on the right hand side 
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Table 1: The positions of the leading Landau singularities generated by the graphs displayed in Fig. llQK in MeV). 5 The 
quantities s a ,si ± denote zeros of the polynomial Pl displayed in Eq. )C.4> . For the mass configuration in Fig. I10l \ 
(Fig.Qot)), s a (s l± ) is given in Eq. OJ (Eq. O})- 











A 






339.6 - Z25.4 


B 






339.6 - Z25.4 


C 






341.6-/26.7 


D 


308.4 


355.8 




E 


308.9 


359.0 




F 


308.9 


359.0 





behaves as follows in the vicinity of the anomalous thresholds, 

CO 

V=yrbi+3y 2 ] =c + 2 c «(^)" /2 > c keC. (C.9) 



«=2 



The polynomial Pi is displayed in Eq. iCAV In other words, the function H generates a singu- 
larity of the square root type at s = si ± (s = s a ) in the case [2] (case [3} listed above, whereas it is 
regular in caseQ] The power of the leading singular behavior oc (Pl) 3 ^ 2 in Eq. (IC91 > agrees with 
the general analysis of these singularities provided in Ref. 13811 . 

It turns out that the topologies displayed in Fig. [6] are the generic ones for the occurrence 
of leading Landau singularities, in the following sense. At two loops, these singularities occur 
iff the four-pion vertex at the inner loop amounts to a charge-exchange one, as is the case in 
Figs. [6j3, C. In Fig.[l0]we display all two-loop graphs that generate leading Landau singularities, 
for both K + — > 3n and Ki — > 3n decays. The solid (dashed) lines denote charged (neutral) pions, 
and the double lines incoming kaons. The graphs in Fig.[T0K-C (D-E) generate singularities in 
the complex plane, at s — s a (on the real axis, at s = si±), see Table[T]for numerical values. We 
conclude that, at two-loop order, leading Landau singularities do occur in all decay channels. 

The results Eqs. (IC.8t . (IC.91 > have important consequences for the decomposition of the am- 



plitudes into regular and singular parts, as proposed in Ref. 112 111 , see Appendix H 



Appendix C.1.3. Singularity at the pseudothreshold 

Finally, we discuss singularities at the pseudothreshold s — s p = {Mk - M3) 2 . As can be 
seen from the explicit expression in Eq. (1B.2I >. lim s ^ J;j y\(s) = -00, whereas y2 stays finite. From 
Eq. dC.8l >. we conclude that H therefore tends to infinity as s — > s p - is. The same is therefore 
true for F(s). 

In summary, we conclude that, in the cases \T\ and [2] the function F(s) is analytic in the low- 
energy region of the complex s-plane, cut along the positive axis. The cut starts at the normal 
threshold s = s t = (M c + M4) 1 and extends beyond the pseudothreshold. In the case [3] the 
cut must be deformed in order to encompass the anomalous threshold at s — s a . We refer the 



reader to Appendix D and Appendix E for more detailed discussions concerning the relativistic 



diagrams. The physical values of F(s) are obtained by approaching the real axis from above. 



5 We use the following masses: M n = 139.6 MeV, M n o = 135.0 MeV, M K = M K+ = 493.7 MeV, M K a = 497.6 MeV. 
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B 



C 




D E F 



Figure 10: The graphs that generate leading Landau singularities in K — > 3rr decays at two-loop order. In Figs. A-C 
(D-F) the singularities are located in the complex plane, at s = s a (on the real axis, at s = *i±), where * denotes the 
invariant mass of the outgoing two rescattered pions, e.g. s = (qi + qi) 1 in graph A. See Tablef^for numerical values. 



Appendix C.2. The discontinuity of F 

The origin of the discontinuity AF (s) = F(s + ie) - F(s - ie) of the function F(s) across the 
cut was identified above, in case of the function H(s): the integration contours for s located on 
the upper and on the lower rims of the cut are not the same. In principle, in the same manner, 
AF(s) could be evaluated from the integral representation Eq. (16.23b . For actual calculations this 
procedure is not very convenient, because T{y, s) is only known in form of the series Eq. (16. 24b . 
For this reason, we now evaluate AF(s) by using a different method. 

Since F(s) and Ai(s) only differ by a constant, see Eq. ( 16.15b . we calculate the discontinuity 
of M{s) in the vicinity of the elastic threshold s = s, and then analytically continue it in s. 
The discontinuity of M{s) in the vicinity of threshold <K A 2 can be obtained directly from 
Eq. ( 16.17b . since only the first denominator is singular: 

. x m C d d l d d k ,, 2 2s N(x) 

x = q 0+ f— + y , - Q + x-^ - +y—-j=-J -M c . 

Calculating the integral over y, we obtain 

AM{s) = - i2 J 2 (<D(g ) ~ $(-<?")) , (Cll) 

where 

<&(y) = a ln( yjb - c-a — y + yjb - c+a - yj — yj(b — c-a — y)(b - c + a — y) 

^ , — , y/ c +( b - c - a -y)+ V c -( fo - c + a - y) fn% ~ 

+ 2flV c + c - m ;= > (C12) 

y/b-y 
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and 

2(M g 2 + M 2 b ) V? (Ml + Mg-MjrQod+g))^ (M a ±M b ) 2 

2M K \Q\ ' 2MdQ| ' C± 2(M2 + M,2)- ( ° 3) 

Closer inspection of the trajectory of y\(s), see Fig. [8] shows that, in order to perform the analytic 
continuation consistently, the variable s in the discontinuity should have an infinitesimal negative 
imaginary part, s — > s - ie. 



Appendix D. Relativistic two-loop integrals - equal pion masses 

The evaluation of the two-loop graphs Fig. [3}\ in the non-relativistic theory is presented 
in Sect. [6] and the holomorphic properties of the pertinent loop functions are discussed in 
Appendix C| The procedure is quite complex and non-standard. We therefore wish to com- 



pare the result with the calculation in relativistic quantum field theory, where the propagators in 
Fig- EH are taken to be the relativistic ones. In this and the following Appendix, we discuss this 
more standard calculation, and, in particular, the holomorphic properties of the pertinent loop 
functions. We stick here to the three topologies shown in Fig. [6] In Appendix D ( Appendix E I, 



we consider the case of equal (unequal) pion masses running in the loops, see Fig.|6A (Figs.[6£, 
C). We use scalar vertices throughout, because derivative couplings do not change the holo- 
morphic properties of the integrals. A comparison between the two approaches is provided in 
|Appendix G| where we show that the non-relativistic calculation approximates the relativistic 
one in a well-defined manner - as it must be for the present framework to make sense. 

Appendix D.l. Notation 

We use dimensional regularization and put 

D 

<o=2~ 2 ' (EU) 
where D denotes the dimension of space-time. Loop integrations are symbolized by a bracket, 

r d D i r d D i r d D k 

<-W ^^^J^- 1 (D - 2) 

We abbreviate Feynman parameter integrals by 

{...h= f dx x {...}, {...}i 2 =2 f dXi f x 2 dx 2 {...}. (D.3) 
Jo Jo Jo 



[^] = §^?(l-^nV (D.4) 



Furthermore, we use the measure 

F(3/2 + co) 4' 
with 

C(W)= (4^' (D ' 5) 
Further, in order to simplify the notation, we will suppress the dependence of the various loop 
functions on the meson masses. 
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Appendix D.2. The loop function K(s) 

As mentioned, we first consider the case where the masses of all pions running in the loops, 
as well as of the outgoing pion with momentum ^3, are taken to be equal [cf. Figs.[3]\ and[6}\], 

M a = M b = M c = M d = M„, q\ = M 2 . (D.6) 

The fact that we consider two neutral pions in one of the vertices does not change the holomor- 
phic properties of the graph, as these pions enter the final expression only through the square of 
their total four momentum, s = {q\ + qi) 2 - The pertinent loop integral is 

G(s\ = [[ - * - )) , (D.7) 



D,D 2 D 3 Z)4 



with 



Dy=M 2 n -k 2 , D 2 =Ml-{Q-kf, D 3 =M 2 n -l 2 , D 4 = M 2 n - (P K - I - k) 2 , (D.8) 
and 

& = (PK-q*f. P\ = Ml, Q 2 = s. (D.9) 
The integral G(s) is proportional to the quantity V121 investigated in great detail for real external 
momenta in Ref . 15711 . In particular, these authors have set up codes to evaluate the ultraviolet 
finite part of G(s) by numerical integration. Here, the purpose is quite different: we wish to 
investigate the holomorphic properties of G(s) as well. Our method to evaluate these integrals is 
therefore necessarily very different from the one developed in Ref. ll57ll . 

We proceed in two steps. First, we identify the part of G(s) that stays finite as D — > 4 
by subtracting the divergent subintegral and by removing the remaining overall divergence. In 
the second step, we express the finite part through a once subtracted dispersion relation in the 
variable s. 

To start with, we consider the inner loop 

J ® = (tTTT ) = C{gS)T{-gS){?\ , z = M 2 - Xl (l- Xl )t, t = (P K -k) 2 . (D.10) 



It is useful to write J(t) in a dispersive manner 05811 . 



/«= I ™ (DM) 



-f 



JAM 2 , °~ ~ 1 

and to single out its divergence at D = 4, 

/(f) = /(0) + /(f), 

such that G(s) in Eq. iD.lb becomes 

G(s)=J(s)J(0) + K(s), K(s)= f — /——J -\ . (D.12) 

JamI o~ \DiD 2 (o- - t) I k 

This integral has still an overall divergence at D = 4. It can be made finite by subtracting its 
value at s = 0, 

K(s) = K(0) + K(s) . (D.13) 

As a result, one has 

GO) = J(s)J(0) + K(O) + K(s) . (D. 14) 

We now discuss the finite part K(s). 
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<1\ 



Figure 11: The triangle graph which corresponds to the amplitude T(s) in Eq. <D . 1 5 1 . We use = (Q + qsY, Q 2 
q\ = Ml,P\ = M\. 



Appendix D.3. Peculiarities ofK(s) 

To illustrate the difficulties one is faced with when investigating the holomorphic properties 
of K(s), we first consider the one-loop integral 7(f) in Eq. ( ID. 101 ). in D dimensions. The ie 
prescription for the square of the masses in the propagators is equivalent to providing the vari- 
able t with a positive imaginary part: the integral Eq. (ID. 101 defines a function J'(t) which is 
holomorphic in the upper half plane Imf > 0, and real on the real axis for t < AM 2 . It may 
therefore be holomorphically continued to the lower half plane through the Schwarz reflection 
principle. The denominator z vanishes at two points X\ t2 6 (0, 1) for t > AM 2 . Therefore, JT(f) 
is holomorphic in the complex plane, cut along the real axis for Re t > AM 2 , and 7(f) is the 
boundary value 7(f) = limi mf _ i o + ST(f). The dispersive representation in Eq. (ID. lit is based on 
these observations. 

This chain of arguments cannot be carried over to K(s) without further ado. To show why 
this is so, we consider the integrand 

rW = {im^l '^ 4 ^'=^-*) 2 CD.15) 

in Eq. (ID . 1 21 > [we drop f in the numerator, because this does not affect the analytic properties of 
T(s)]. The quantity T(s) corresponds to the triangle diagram displayed in Fig. QT] 

The integral is ultraviolet convergent and can thus be evaluated at D — 4. After integration 
over the momentum k, we are left with the Feynman parameter integral 

T(s) = [z~% , Z = x\Ml + cr(l - x 2 ) - x 2 x x [s{\ - xi)x 2 + (1 - x 2 )(M 2 K - M 2 K )} - ie . (D.16) 

We have explicitly displayed the ie prescription. In this case, we cannot replace ie by providing s 
with a positive imaginary part: the denominator z vanishes at two points on the xi — 1 boundary, 
for any value of s, provided that Mk > Vo 1 + M„, a condition which is met for the physical 
value of the kaon mass. In other words, at e = 0, the integral Eq. ( ID. 16b is not well defined for 
real, physical values of M 2 ., M 2 K , at any s, and cannot straightforwardly be considered to be the 
boundary value of a holomorphic function T(s). This renders the discussion more complicated 
than in the case of J(s) M59I1 . A similar remark applies to K(s). 

A way out of the problem is based on the following observation. First, examining Eq. ( ID- 1 6b . 
one notes that T(s) has a unique extension T(s) into the product of the upper half co mplex p lanes 
of s and M 2 K One then considers the kaon mass to be an adjustable parameter [45, 59, 60|]. For 



6 It was already noted in Ref. |39i] that T(s) is holomorphic when all three variables M 2 K , s, and q 2 lie in the same half 
planes. 
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sufficiently small value of M 2 K , the above-mentioned problem does not exist - the function 7~(s) 
satisfies an unsubtracted dispersion relation in the variable s, with threshold s, = AM 2 . For 
physical values of the kaon mass, T(s) is then defined by analytic continuation of the dispersive 



representation through the upper complex M^.-plane, and T(s) = ]imj ms - t o+,imM K ^o+ T(s) [45 



59L160D- 

This procedure may now be directly carried over to K(s) - the final integration over <x in 

Eq. ( 1D.121 I does not affect the arguments just outlined, aside from the necessity to use a subtracted 

1 

A'' 



dispersion relation. Concerning the continuation in M 2 we consider the Cauchy representation 



■«n Jam 1 - x ( x ~ s > 

where C denotes the path of integration in the complex s-plane, and Ai<C stands for the dis- 
continuity of i{ across C. For small kaon masses, C may be taken to run along the real axis 
from AM 1 , to infinity, and 7C is holomorphic in the complex s-plane, cut along C. The discon- 
tinuity is itself a holomorphic function t^K(s) in a certain region of the complex s-plane, with 
singularities Si(Mk, M„) whose positions depend on the pion and kaon masses. These singular- 
ities generate singularities in the function 'K on higher Riemann sheets. During the process of 
increasing the value of the kaon mass to its physical mass, one has to make sure that the sin- 
gularities s, do not cross the path C, by eventually deforming it properly 11451 159[ |60[ - in other 
words, the singularities in 'K may eventually move to the first Riemann sheet. Once the contin- 
uation has been accomplished, it is straightforward to read off the analytic properties of "7C(s), 
together with e.g. its threshold behavior, which is of particular interest in our case. Finally, one 
has K(s) = lim Ims ^ 0+ , ImM ^o+ 'K(s). 

Appendix DA. The discontinuity A'K 

It remains to evaluate the discontinuity in the dispersive representation Eq. ( ID. 171 >. For suf- 
ficiently small values of the kaon mass, unitarity of the S -matrix fixes it unambiguously. It is 
given by the angular integral 

P {z) = ( ! - ZTfii) ' A ^\[ M i + 3M ^- s ] ' B^^-A l l 2 {s,M 2 K ,Ml). (D.18) 

After the substitution (J3(z) - l)/(J3(z) + 1) = u, the integral can straightforwardly be performed. 
For the following discussions, it is useful to introduce the functions 

0(x,y;z + ,z_;m) = m 2 (r 2 + - r 2 ) - {x+y)z+r+ + (x -y)z_r_ , r ± = ln(l + z±) - ln(z ± - 1) , 
F(x,y;z+,Z-;m) = €>(x,y;z+,z_; m)\ r± ^ R± , R ± = ln(l +z±) - ln(l - z ± ) , 

G(x, y; z + ,Z-;m) = 2m 2 (R- +R + )-(x + y)z + -(x- y) Z - . (D. 1 9) 



7 Here and below, we denote the independent variable in the discontinuity with s. 



38 



The discontinuity for small values of the kaon mass is given by 

AK(s) = ^£s[6 + B- l D(s)] , D(s) = <t>(A,B;W + ,W-;M n ), W ± = |l - i^J . 

(D.20) 

According to the discussion above, this representation must be continued analytically to the 
physical value of M 2 K . Potential singularities in the discontinuity occur for those values of s 
where the arguments of W± vanish, or B vanishes. To investigate the first possibility, we observe 
that the quantity (A - AM 2 ) 2 - B 2 agrees up to a factor with the polynomial B^ in Eq. ( IF. 31 ) that 
determines the location of potential leading Landau singularities in the graph Fig. [6}\, 

AM 2 1 
(A-4M 2 ) 2 -fi 2 = — ^(s-s,) 2 , s, = -(M 2 -M 2 ). (D.21) 

The zeros of B occur at s - AM 2 , (Mk + M„) 2 . Because M 2 K is equipped with a positive imaginary 
part, the zeros at s = su (Mk ± M„) 2 are located above the path C of integration. At s = AM 2 , 
the discontinuity vanishes. 

There are thus four intervals to be distinguished, 

Interval Illb Ilia II I 

[s,,s{] [Si,S p ] [s p ,s c ] [s c ,oo) 

where 

s f =4M 2 , s p = {M K -M n f, s c = {M K +M n f. (D.23) 

We now list the discontinuity in the four intervals for physical values of the kaon mass, and use 
the notation 

v ( 5 ),w . f 2 w2x,l/2 



B=^\A{s,M' K ,Mt t )Y l \ (D.24) 



Interval I [A7C is imaginary] 



Here, B — B, and W± - 1 > 0. The discontinuity may therefore simply be evaluated from 
the expression given in Eq. {D.20) . As a check we have verified that the expression for A'K, 
evaluated at Mk = M„ = 1, agrees with the discontinuity of the two-loop function V\ worked out 
in Ref. JH, Appendix A.2]. 

Interval II [A'K is imaginary] 

The function D(s) is evaluated as in interval I, with B - -iB. 
Interval HI [A'K is complex] 

In this interval, B = —B, and the discontinuity is obtained from 

1/2 

-4 IV! ~ 

DO) = inG(A, B; W+, W-, M„) + F(A, B; W+, W-,M n ) , W ± =c ± 



1 - 



A + B 



(D.25) 



The numbers c ± are 

Interval Ilia Illb 



c + 1-1 (D.26) 

c_ 1 1 
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0.30 0.35 0.40 0.45 0.50 



Vs [GeV] 

Figure 12: The discontinuity A'K'(s) as a function of t/s, in the equal mass case, in arbitrary units. The solid (dashed) 
curve stands for the imaginary (real) part of A'K. The singularity at the pseudothreshold -^sj, is clearly seen. 

Ims 

S t Si s p 
L--., J -, • - r • Re s 

c 



Figure 13: The path C in the dispersive integral Eq. (D. 17K for physical kaon mass. The symbols s,, Si, and s p are 
defined in Eq. jD.23t . 

The discontinuity develops a singularity at the pseudothreshold s = s p . This is due to the fact 
that the angular integral in Eq. (ID. 18) must be deformed to infinity in the complex plane as s 
approaches s p J45ll . and the angular integral finally diverges at s — s p . On the other hand, there 
is no singularity at s — s\ [second order zero in Eq. (ID.21H . 

In Fig. [12] we display the discontinuity A'K(s) as a function of Vs. The solid (dashed) 
curve stands for the imaginary (real) part of MK. The singularity at the pseudothreshold s p is 
clearly visible. From the explicit expression of the discontinuity, one furthermore concludes that 
it generates a square root singularity at threshold, 

AK(s) = v(s)d(s) , (D.27) 

with d(s) holomorphic at s = 4M^. 

Appendix D.5. Holomorphic properties ofK(s) 

We deform the path C in the representation Eq. (ID. 17t such that the singularity of the discon- 
tinuity at the pseudothreshold s p is avoided in the manner indicated in Fig. [13] as is requested by 
the continuation in M 2 K . The function "Kis) is thus holomorphic in the complex s-plane, cut along 
the real axis for s > AM 2 . Furthermore, on the upper rim of the cut, it is holomorphic as well. 
The singularity of the discontinuity at s — s p shows up in 'K(s) only on the lower rim of the cut, 
where it diverges when s p is approached from below. We may perform the dispersion integral 
Eq. (ID. 171 > numerically. The result is displayed in Fig. [14] with the variable s at the upper rim of 
the cut. 
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0.3 



Figure 14: The function 'K(s) evaluated numerically from Eq. JD.17I . in arbitrary units, with s at the upper rim of the 
cut. The solid (dashed) line denotes the real (imaginary) part of 'K(s). 



Appendix E. Relativistic two-loop integrals - unequal pion masses 

We now investigate the graphs Fig. |6j3, C, see also Fig. |3]\. The corresponding ultraviolet 
finite loop functions 'Kb.c are again written as a dispersive integral, 



KM 



2ni 



C f*oo 
J(M r H 



+M d ) 2 x ( x ~ s ) 



I = B,C- 



(E.l) 



The discontinuities are given by angular integrals 



Wis) = 



ivcd(s) 
256tt 3 J_ 



BAz^n ^l \ +l\dz; s e [(M c + M d )\ co| . 
Biiz) + 1 



where 



Vcd(s) = 



A ll2 (s,M 2 c ,M 2 ,) 



Bb(z) = 1 - 



Ml 



A+Bz 



1/2 



Bciz) = 



M 



2 \ 



1 - 



A + Bz) 



1/2 



Vcd(s) 



A = M z K + Mf,--(s + Mt-M z d )(M' K + s-qi), B 



A 1/2 (s,M 2 K ,qj). (E.2) 



For later use, we introduce 



B 



Vcd(s) 

2s 



UO,M|,^)| 



2 1/2 



(E.3) 



Appendix E.l. Case I 

We first consider the graph in Fig. |6ji, 



M a = M h = M c = M n , M d = Mjs , ql^M 2 , 



(E.4) 



Singularities are generated by the zeros in 



W + 



4M. 
A + B 



2 X 1/2 



(E.5) 
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The quantity (A - AM 2 ) 2 - B 2 agrees up to a factor with the polynomial in Eq. (1F.31 >. which 
shows that corresponding singularities in A'Kb occur at 

*i± = \ [Ml + 2M 2 , - 3M 2 ] ± ^- [(M 2 - M(M|, M\, AM 2 n )\ 12 . (E.6) 

We have to distinguish the following intervals: 



Interval IIIc Illb Ilia II I 

[s t ,si-] [s l+ ,s p ] [s p ,s c ] [s c ,oo) 



where 



s t = (M„ + M n f , s p = (M K - M„of , s c = (M K + M,fi) 2 . 
The discontinuity is 



where £>b(s) is given by the following expressions in the intervals I— III: 

Interval I [A'Kb is imaginary] 
Here, one has B = B, and 



Db(s) — <&(A, B; W+, W-\ M n ) , W ± = \l 



4M- 
A + B 



2 vl/2 



(E.7) 



(E.8) 



(E.9) 



(E.10) 



Interval II [A'Kb is imaginary] 

The discontinuity is evaluated as in interval I, with B = -iB. 

Interval HI [A'Kb is complex] 

In this interval, one uses B = —B, and 

D B (s) = inG(A, B; W+, W-,M n ) + F(A, B\ W+, W-, M„) , W ± = c ± 
where 

Interval Ilia Illb IIIc 
c + 1 —i -1 

c_ 111 



4Ml 



A + B 



1/2 



(Ell) 



(E.12) 



Appendix E.2. Holomorphic properties ofKs 

The discontinuity A'Kb develops a singularity of the square root type at s — s\ ± . The path of 
integration in the dispersive representation Eq. (IE. 1) is indicated in Fig. [15] The function 'Kb is 
holomorphic in the complex s-plane, cut along the positive real axis for s > (M„o + M^) 2 . It is as 
well holomorphic at the upper rim of the cut, and diverges at the pseudothreshold s = s p when 
approached from below the cut. Further, at s — s\±, it has singularities of the square root type on 
the lower rim of the cut. Its precise behavior there is indicated in Eq. ( 1C.9I ). 
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L---, • • n * r Re S 



Figure 15: The path C in the dispersive integral Eq. dE. It . with mass assignments as displayed in Eq. jE.4K for physical 
kaon mass. The symbols s t , si + , and s„ denote threshold, two anomalous thresholds, and the pseudothreshold, 
respectively, and are given explicitly in Eqs. JE.61 . jE.81 . 



Im ,S 



4M^ 



c 



Rc s 



Figure 16: The singularities s„, s* in the discontinuity A'Kc(^), for the mass assignments Eq. iE. 131 . as a function of the 
kaon mass (dotted lines). The symbol s a is defined in Eq. <F.4t . .s* denotes its complex conjugate, and the arrows indicate 
the direction of increasing kaon masses. The singularity s a intrudes into the original path of integration, such that C must 
be deformed [dashed line], see also Fig. 1171 



Appendix E.3. Case II 

Finally, we consider the case where the pions in the inner loop are neutral and all the others 
charged, 

M a =M b =M n o, M c — Mj — M K , q\=M 2 n , (E.13) 

see Figs. [3}\, |6p. This case is the most intriguing one. Performing the angular integral in 
Eq. (lE.2t and investigating the singularities as before, we find that the ones in the pertinent 
square roots W ± are given by the polynomial ijgj, displayed in Eq. (IF. 3 1 > . We display the lo- 
cations of the singularities s a , s* a as a function of the kaon mass in Fig. [16] [dotted lines]. We 
see that, as the kaon mass is increased, the singularity s - s a intrudes into the original path of 
integration. As a result, the path C must be modified [dashed line], see also Fig. [17] 



Im s 



Res 



' • ' S a 



Figure 17: The path C in the dispersive integral Eq. JE.lt , with mass assignments as displayed in Eq. IE. 13t . for physical 
kaon mass. The symbols s t , s a , and s p are defined in Eqs. ID. 231 and )F.4t . 
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Appendix E.4. Holomorphic properties ofKc 

From the above discussion, we conclude that 'Kc(s) is holomorphic in the complex s-plane, 
cut along the path C displayed in Fig. [T7] It is holomorphic at the upper rim of the cut, and 
diverges at the pseudothreshold s p when approached from below. Further, at s — s a , it develops 
a singularity of the square root type on the lower rim of the cut. Its precise behavior there is the 
one indicated in Eq. dC.9l l. The singularity at s* is located on the second Riemann sheet. 



Appendix F. Leading Landau singularities 

Some of the graphs that correspond to Fig. [3]\ have the property that they generate singu- 
larities that correspond to the solutions of the pertinent leading-order Landau equations. These 
singularities occur because the kaon is unstable. The positions of these singularities are given by 



the two solutions of a polynomial equation of second order in s [57] 



P L (s, q\, M\; M a , M b , M c , M d ) = , (F. 1) 



with 



P L = A(s, M 2 , M 2 c )A(ql (M a + M b ) 2 , M 2 ) - Q 2 (s, q\, M 2 K ; M a , M b , M c , M d ) , 
Q = 2M 2 (M 2 K -s-q\) + (qj - (M a + M b f + M 2 )(s + M 2 - M 2 .) . (F.2) 

For the graphs Figs|6]\-C, this polynomial becomes 

P L (s,M 2 n ,M 2 K ;M n ,M„,M„,M n ) = -16M* [s - si] 2 , 
^ = P L ( 5 ,M 2 ,M 2 ;M„M ff ,M„M,o) = -16M 2 M^ [s - s l+ ] [s - Sl _] , 

= P L (s, M 2 -, M\; M^,M^,M n , M„) = -16M 2 M 2 [ s - s a ] [s - s* a ] . (F.3) 

Here, s* a denotes the complex conjugate of s a . The thresholds are 



Sl = ^(M 2 K -M 2 n ), 



si ± = \ [M 2 k + 2M 2 - 3M 2 ] ± jL- [(M 2 n - Ml )A(M 2 K ,MlAMl)] 112 , 
*a = \ K + 3M ' - AMl A ~ 2M^ V M l - M ^M\, Ml, AM 2 *)] 1/2 . (F.4) 

The positions of the leading Landau singularities of the triangle graph displayed in Fig. [18] are 
obtained from a polynomial equation Pj(s,q\, M 2 K ; Mx, M c , Md) = Il37ll . It is interesting to 
note that, for Mx - M a + M/„ the polynomials Pi and Pj are proportional to each other, and 
the positions of the leading Landau singularities in the two-loop graphs Fig.[3}\ and the triangle 
graph Fig.[l8]thus coincide. This corroborates the statement made in Ref. 1 59, Section 3] on the 
origin of the singularities in the graph |6]\. 

The discussions in the previous appendices show that the discontinuities of the graphs in 
Fig- Si, C do have a singularity at *i±, s a , s* a , while the one of Fig.|6}\ is holomorphic at S\. As 
a result of this, the pertinent loop functions generate square root singularities at s\ ± , s a on the 
lower rim of the cut. 
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<1\ 



Figure 18: The triangle graph with general masses. The positions of its leading Landau singularities agree with those of 
the two-loop graphs in Fig.|5]\ for M x = M a + Mb, for general q\. The notation is Pt = (Q + qsf, Q 2 = s, P\ = M\. 



Appendix G. Relativistic vs. non-relativistic approach 



As already mentioned in Appendix D the non-relativistic calculations should reproduce the 
relativistic results graph by graph in a systematic expansion in e - otherwise, the whole frame- 
work does not make sense. In other words, it should be possible to adjust the couplings in the 
non-relativistic Lagrangian at order e" so that the relativistic expression of any graph, expanded 
in powers of e, up to and including 0{e") would be reproduced by a uniquely identifiable set of 
non-relativistic graphs. 

In this Appendix, we check this equivalence for the relativistic two-loop graph, considered 
in | Appendix D| To this end, we first compare the expression for the non-relativistic discontinu- 
ity, Eq. ( 1C.1U . to the expression obtained in the relativistic theory. For simplicity, we restrict 
ourselves to the case M a - Mb = M c - Mj = M3 = M„. The comparison is carried out in 
the threshold region Illb defined in Eq. (1D.22I ) where the non-relativistic expression has been 
originally derived. The discontinuity in the whole non-relativistic region is obtained by using 
analytic continuation. 

The quantities A and B, which were defined in the relativistic theory, see Eqs. (ID. 1 81 > and 
( ID. 241 ). are related to the quantities a, b, qo through 

4M£ A B 2M K \Q\ 

a- , b-—, qo = — , k- — — . 

k k k 

From the above equation we find (cf. Eq. (1D.20I >)) 

»'/ 2 / 4 Ml v 2 

and 



/ a \ [ / 2 / 4M 2 n \'/ 2 _ 

1-- = 1 ^ =W ± (G.2) 

I b + qo) \ A + B) 



0)( ?0 ) - ®(-<Zo) = TTj^: (2Ml(R- - R + ) + (A + B)W + - (A - B)W.) , 

/? ± = ln(l + W ± )-ln(l - W ± ). (G.3) 

The final expression for the non-relativistic discontinuity is obtained by substituting Eq. ( lG.3t 
into Eq. ( IC.lll l. It is seen that this expression differs from the pertinent relativistic expression 
given by Eqs. dD~20l i. (EH), dD~25l ) by 

Disc,(i)|^ - EHsctoU = ^)P(s) , P(s) = ^3 (6 - |) , (G.4) 
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B 



Figure 19: Non-relativistic diagrams that remove the right-hand side of Eq. <G.4t . The thick squares stand for the 
counterterms in the non-relativistic K — > 3rr Lagrangian at order a. 



where v(s) = (1 — 4M^/i)^ 2 , and F is defined in Eq. (ID. 19b . One may ensure that the quantity 
P(s) is a low-energy polynomial with real coefficients. Consequently, the relativistic and non- 
relativistic amplitudes themselves differ by iv(s)P(s) + P'(s), where P'(s) is another low-energy 
polynomial. 

In order to establish the origin of the above polynomials, let us first compare the relativistic 
and non-relativistic one-loop integrals in the vicinity of the threshold. As seen from Eq. (12.91 l. 
only the imaginary part of the one-loop integral above threshold survives in the non-relativistic 
approach (the expression below threshold is obtained through analytic continuation). This cor- 
responds to the following replacement in the expression of the one-loop integral J(s) above 
threshold 

v(s) - 1 1 - v(s) 

v(s) In — + 2 = vO) In — + 2 + wv({) -> inv(s) . (G.5) 

v(s) + 1 1 + v(s) 

The dropped part is a low-energy polynomial with real coefficients (in the whole non-relativistic 
region) and corresponds to adding the counterterms to the K — > 3n vertex shown in Fig. [19R. 
The numerical values of these counterterms are fixed at order a. 

It can be straightforwardly checked that inserting the vertices with the counterterms inside the 
pion loop, as shown in Fig.[T9B. and evaluating this diagram in the non-relativistic theory yields 
exactly the result iv(s)P(s). Consequently, this term takes into account the renormalization in the 
non-relativistic theory at 0(a). Finally, the polynomial P'(s) corresponds to the renormalization 
of the K — > 3n vertices at 0(a 2 ). At the order we are working, only tree diagrams with these 
counterterms contribute. 

To summarize, we have demonstrated that all differences between the relativistic and non- 
relativistic amplitudes in the low-energy region can be removed by changing the renormalization 
prescription in the K — > 3n vertices. Consequently, relativistic and non-relativistic theories are 
physically equivalent at this order. Obviously, this also holds for the expansion of the pertinent 
Feynman graphs in powers of e. 



Appendix H. The decomposition M = Mo + v cd M\ 



One of the crucial ingredients in the work of Cabibbo and Isidori 112 111 is the decomposition 
of the decay amplitudes M into the form 

M = M + v cd (s)M l . (H.l) 

The amplitudes Alo.i are assumed to be analytic in the physical region of the decays, with square 
root singularities at the border of the Dalitz plot, associated with different nn thresholds ||2|]. 



46 



In this Appendix, we comment on this decomposition, and start the discussion with a simpler 
case, the relativistic one-loop integral 



M 2 -I 2 M 2 d -(l- P) 2 



■/„/<*) = ( ,/ „ ,„ ;, , s = P 2 . (H.2) 



Let 



Jcd(s) = Jcd(s) - Jcd(O) . (H.3) 

The function J c d is analytic in the complex s-plane, cut along the real axis for s > (M c + Md) 2 . 
Its explicit form near the threshold s = (M c + Md) 2 is indeed of the form Eq. ( IH.ll i. Define Mo 
through 

Us) = Mo + . (H.4) 

Ion 

It then follows that Mo has no right hand cut, because the term i v c d(s)/\6n reproduces the 
discontinuity of J(s) at s > (M c + Md) 2 - On the other hand, this term develops singularities at the 
pseudothreshold s = (M c - Md) 2 , and at s - 0. Because J c d is analytic there, the amplitude Mo 
develops the same singularities, with opposite sign, such that the sum is regular at these points. 
Since s = 0, s — (M c - Md) 2 are outside the kinematic region of interest in the present case, these 
singularities do not matter, and one concludes that J c d(s) indeed has a decomposition Eq. dH.U 
with amplitudes Mo,i that enjoy the properties proposed in Ref. H. 

For the two-loop graphs discussed above, the situation is more complex. Let us first discuss 
the situation in case of the function H(s) introduced in Eq. (IC.8b . with mass assignments that 
correspond to the diagram Fig. |6ji. We start from the discontinuity 

AH(s) = H + (s) - H-(s) = — — -vFyiCyi + 3y 2 ) , s>{M„ + M 7fi ) 2 , (H.5) 

worked out in Eq. (IC.8I 1. We note that y\ = -v 2 fi(v 2 ),y2 = — /2(v+ ) near threshold, with f,(z) 
holomorphic at z = 0, and /)(0) > 0. Now define the amplitude Hq(s) through 



1 



H(s) = -AH(s) + H (s) , s > (M nB + M n f . (H.6) 

Evaluating the discontinuity on both sides, it is seen that Hq is holomorphic near threshold. 
Therefore, 

H(s) = Hois) + v+oHds) , Hi = —AH(s) , (H.7) 

2v +0 



with Ho.i holomorphic at threshold. However, as is shown in Appendix C.1.2 AH is singular 
at s — s\ ± . Because His) is regular there, Ho must be singular as well. A completely anal- 
ogous argument may be used to discuss the inadequacy of the decomposition ( 1H.U in case of 
the relativistic integral discussed in Appendix E. 1 Numerical values for the positions of these 



singularities are provided in Table [T 

We find it instructive to display the singular behavior of Mo, 1 in a fully explicit manner, in a 



simplified example inspired by the relativistic integrals discussed in Appendix E Consider the 
dispersive integral 

1 c Cdx^x-AMlyfsT- 
A J " 



/(*) = - 2 n —^- , siC. (H.8) 
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Im x 



4M= si_ Sl _+A 
L , I F J ► Re :r 

c 

Figure 20: The path C along which the integral is taken in Eq. dH.8> . 



Here, -<Jx - 4M 2 stands for the standard phase space factor, - x is a substitute for the 

square-root singularity in the discontinuity at s — s\-, and A denotes a (mass-scale) 2 specified 
below. The path C is indicated in Fig.[20l compare with Fig. [151 

The function I(s) is holomorphic in the complex s-plane, cut along the real axis for s e 
[4M 2 , s\- + A]. Further, 7(s) can be holomorphically continued through the upper rim of the cut, 
because the singularity at s = si_ is not effective there - the path C can be suitably deformed. On 
the other hand, approaching the cut from below, it is seen that C cannot be deformed any more at 
s = si_, as a result of which one encounters a (pinch) singularity there. 

To perform the integration, we set A = s\- - 4M 2 , x = 4M 2 + Ay, 



m I 

Jo 



dy Jy Jl-y + ie s - 4M 2 
I( S ) = G(z) = lim V V , z = t— S . (H.9) 

y-z A 



We find 

r l dy JyJl-y f 2 dyJyJy-l 
G(Z)= I _ +ij . _ =Gi<: ■)+/(■;:(;), (11.10) 



with 



4VzVl -z 



Gi = zVr ViVl - z + pi(z) , Gi - Vz V 1 - z arctan ^+p 2 (z), 

(3z - 2) V2 

pi=-a(z-i), p 2 = V2 + ln(3+2V2)(z-i), (H.ll) 

where z is located on the upper rim of the cut, near z — 0. One concludes that 

/(i) = i ^s-AMlh{s) + h(s) , A = | ^/sT~^ , 7 = pi(z) + iG 2 (z) , (H. 12) 

with 7o,i holomorphic at threshold. This is indeed a decomposition of the type Eq. ( IH.ll i. Note, 
however, that 7i develops a square root singularity at s = s\— As was said above, I(s) is holo- 
morphic on the upper rim of the cut. Therefore, 7q develops a singularity at s = si_ as well, 
canceling the one generated by 7i . Indeed, near z = 1 , one has 



g 2 = VzVT^I 



4 VzVl -z 
arctan — - n 

(3z - 2) V2 



7>2(z), (H.13) 



as a result of which the singularities in 7o,i cancel in the sum ( 1H.12I ). On the other hand, on the 
lower rim of the cut, near s - s\-, the singularities add, and I(s) shows a square root behavior. 
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